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Aims. Archeops is a balloon-borne experiment widely inspired by the Planck satellite and by its High Frequency Instrument (HFI). It is mainly 
dedicated to measure the Cosmic Microwave Background (CMB) temperature anisotropies at high angular resolution (~ 12 arcminutes) over 



1^ a large fraction of the sky (around 30 %) in the millimetre and submillimetre range at 143, 217, 353 and 545 GHz. Further, the Archeops 

O 353 GHz channel consists of three pairs of polarized sensitive bolometers designed to detect the polarized diffuse emission of Galactic dust. 

^ Methods. We present in this paper the update of the instrumental setup as well as the inflight performance for the last Archeops flight campaign 

CZ3 in February 2002 from Kiruna (Sweden). We also describe the processing and analysis of the Archeops time ordered data for that campaign 



which lead to the measurement of the CMB anisotropies power spectrum in the multipole range ^ = 10 - 700 (Benoit et al. 2003a, Tristram 
et al. 2005) and to the first measurement of the dust polarized emission at large angular scales and its polarized power spectra in the multipole 
range C = 3-70 (Benoit et al. 2004, Ponthieu et al. 2005). 

Results. We present maps of 30 % of the sky of the Galactic emission, including the Galactic plane, in the four Archeops channels at 143, 217, 



. 353 and 545 GHz and maps of the CMB anisotropies at 143 and 217 GHz. These are the first ever available sub-degree resolution maps in the 
- - ' millimetre and submillimetre range of the large angular-scales Galactic dust diffuse emission and CMB temperature anisotropies respectively. 
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1. Introduction 

The measurement of the Cosmic Microwave Background 
(CMB) anisotropies in temperature and polarization is a funda- 
mental proof of modern cosmology and of the early Universe 
physics. Since the first detection of the CMB anisotropies 
by the COBE satellite in 1992 ( |Smootetal. (1992)t , a 
large number of ground-based and balloon-borne ex- 
periments such as DASI (pialverso n et al. (2002)|, CB I 
(Mason et al. (2003)|, VSA ([Dickinson et al. (2004)^ 
BOOMERanG (|Netterfield et al. (2002)t or Maxima 
( |Hanany et al. (2000)> have measured the CMB angular 
power spectrum from a few-degrees down to sub-degree 
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scales. However, simultaneous observation of very large 
and small angular scales have proved to be particularly 
difficult, as it requires both large sky coverage and high 
angular resolution. This has been achieved, first, by Archeops 
( [Benoit et al. (2003a)| [Tristram et al. (2005b)| l which has 
measured the CMB power spectrum in the multipole 
range 10 < ^ < 700. Then, the WMAP satellite mission 
( [Bennett et al. (20(n)| ) has detected the CMB anisotropies, 
both in temperature and polarization. 

Archeops, described in details in [Benoit et al. (2002)[ is a 
balloon borne-experiment designed as a prototype for the High 
Frequency Instrument (HFI) of the Planck satellite. Its tele- 
scope and focal plane optics are widely inspired by the Planck 
design. The implementation of the measurement chains: cryo- 
genics, optics, bolometers, readout electronics was a successful 
validation of the innovative design. Further the data processing 
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was a learning process for future members of the HFI team. 
Archeops performs circular scans on the sky with its optical 
axis tilted 41 degrees with respect to the horizon by spinning 
the gondola at 2 rpm. This scanning, that is combined with the 
proper motion of the Earth, leads to 30% sky coverage in about 
12 hours of flight. With a rotating gondola, the Sun above the 
horizon produces a dominant parasitic signal. The optimal way 
of avoiding this, while having the longest integration time, is 
by having a long duration (Arctic) night-time balloon flight. 

The Archeops payload was successfully launched 
three times. First, from Trapani (Italy) in July 1999 
l |Benoit et al. (2002)) for a 4 hours test flight. Then from 
the Swedish Esrange station (near Kiruna at 68 deg. latitude 
North, just above the Arctic circle) operated by the French 
Centre National d'Etudes Spatiales (CNES) and the Swedish 
Space Corporation in January 2001 (hereafter KSl flight) 
for a 12 hours flight, and finally in February 2002 (hereafter 
KS3 flight) for a 24 hours flight from which 12 night-hours 
were exploited for scientific purposes. In the KS3 flight, 
a stratospheric altitude of 35 km was reached, reducing 
significantly the contamination from atmospheric (mainly 
ozone) emission with respect to ground-based measurements. 
Additional information about the Archeops flights may be 
found at our web-site'. 

The Archeops bolometers are grouped in four frequency 
bands at 143 GHz (8 bolometers), 217 GHz (8 bolometers and 
a blind one), 353 GHz (3 polarized bolometer pairs), 545 GHz 
(1 bolometer). The 143 and 217 GHz channels are dedicated to 
the measurement of the temperature angular power spectrum of 
file CMB (Benoit et al. (2003a)| [Tristram et al. (2005b)> . The 
353 and 545 GHz channels allow the monitoring of both at- 
mospheric emission and Galactic thermal dust emission. In ad- 
dition, the polarization of the diffuse Galactic dust emission has 
been measured for the first time using the 353 GHz polarized 
bolometers jBenoit et al. (2004)||Ponthieu et al. 92005)> . 

We present here the processing of the Archeops data 
for the KS3 flight going from raw time ordered data to 
maps of the sky. The Archeops data processing was specif- 
ically designed to cope with the characteristics of the 
scanning strategy and has similarities with Planck-HFI data 
processing. Other dedicated processing techniques are de- 
scribed in details in Hanany et al. (2000)' "Lee et al. 20011 
Stompor et al. 2002| Rabii et al _2005 Ruhl et al. 20031 
Masi et al. (2005)| IKuo et al. 2004) to deal wifli the Maxima, 
BOOMERanG, ACBAR spider web bolometer-experiments 
data and with WMAP HEMT all sky survey satellite 
dati jHinshaw et al. (2003)| 

This paper is organized as follows. Sections |2l and |3l de- 
scribe the instrumental set-up and the in-flight performances 
of Archeops during the KS3 flight. Section|4]presents the pre- 
processing of the Archeops data. In section [S] we describe the 
offline pointing reconstruction. Section|6ldeals with the optical 
and time responses of the instrument. Sections 0and|8lpresent 
the characterization and treatment of systematics and noise in 
the data. In section |5] we discuss the intercalibration and ab- 
solute calibration of the Archeops data. Finally, section [Tol 

' [http : //www ■ archeops . org| 



Fig. 1. Scheme of the Archeops gondola 



presents the construction of the Galactic and CMB Archeops 
maps. We conclude in section [TTI 

2. Technical description of the experiment 

In this section we describe the main aspects of the instrumental 
setup of the Archeops experiment. Particular interest is paid to 
changes performed on this setup since the Trapani test flight 
( IBenoit et al. (2002)> . 

2.1. Gondola 

The Archeops gondola was redesigned since the Trapani flight, 
in order to gain some weight and try to reduce the main para- 
sitic effect that was observed then. The most likely interpre- 
tation was indeed that stray light reflecting from or emitted 
by inhomogeneities on the balloon surface was the main cul- 
prit for this large scale parasitic effect. The new baffle contains 
lightweight highly reflective material (Fig.^ in a staircase-like 
layout so that the entrance of the gondola is highly "reflective" 
for downward rays. To further reduce the systematics effects, 
between flight KSl and KS3, the engine driving the gondola 
spinning was moved from the top of the gondola (and rigidly 
fixed to it) to an higher location along the flight chain, 60 meter 
above the gondola. This allowed to strongly reduce non station- 
ary noise induced by the swivel engine. 

2.2. Attitude control system 

The Archeops attitude control systems are composed of gy- 
roscopes, a GPS and a Fast Stellar Sensor. The gyroscopes 
are unchanged since Trapani flight. A high precision z-axis 
laser gyroscope based on the Sagnac effect was added for long 
term relative azimuth reconstruction. This is needed for day- 
light data, when detected stars are not enough to track rotation 
speed changes. The GPS was changed since Trapani because 
of a failure at high altitude. A 1-m diameter circular loop of 
copper wire was added and used as an Earth magnetic field de- 
tector to perform a rough (5 degree accuracy) absolute azimuth 
reconstruction. The electromagnetic influence of the pivot rotor 
and its associated wire on the loop signal disappeared from the 
KSl to the KS3 flight, as the rotor was moved upwards along 
the flight chain. 

The Fast Stellar Sensor (ESS) is a 40-cm optical telescope 
equipped with 46 photodiodes mounted on the bore-sight of 
the primary mirror (Fig. Q for a-posteriori accurate point- 
ing reconstruction. The photodiodes are aligned along a line 
which is perpendicular to the scanning direction. Each photo- 
diode covers a 7.6 arcminutes (parallel to the scanning direc- 
tion - para scan) by 1.9 arcminutes (perpendicular to the scan- 
ning direction - cross scan) area on the sky. The ESS sweeps 
a 1 .4 degree wide ring at constant elevation during a payload 
revolution and its center is mechanically within one degree of 
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Fig. 2. Optical layout of the Archeops focal plane. 



the main submillimetre telescope pointing direction. The FSS 
has been improved with respect to the Trapani configuration 
l |Benoit et al. (2002)) : a red filter has significantly diminished 
the background and parasitic noises have been suppressed. A 
fuU report on the FSS is given by psfati et al. (2003)) During the 
flights, about 100 to 200 stars per revolution were detected by 
the FSS. A detailed description of the pointing reconstruction 
can be found in Sect.|5] 

2.3. Detectors 

For the last KS3 flight campaign, the detection was ensured by 
an array of 21 spider web bolometers (B ock et al. (1996)| l of 
the same type as for the Maxima (100 mK) and BOOMERanG 
(300 mK) experiments. For each bolometer, a Neutron 
Transmutation Doped Germanium thermistor is fixed on a sili- 
con nitride micromesh designed to absorb submillimetre light. 
The bolometers are cooled down below 100 mKby an ^He/"^He 
dilution cryostat ( Benoit & Pujol (1994)j and were optimized 
for the expected background loads at this temperature, varying 
from 2 to 8 pW depending on frequency. The bolometers were 
built at JPL/Caltech in the context of the development of the 
Planck HFI instrument ( |Lamarre, J.M. et al. (20 03) I. 

Bolometer characteristics were measured from standard I- 
V curves obtained during ground-based calibrations at zero 
power load, after the bolometers were made blind (see sec- 
tion I3.2.2I I. In order to prevent radio frequencies contamina- 
tions, each bolometer is kept in a copper A/4 cavity acting as a 
Faraday cage for maximal absorption. 

The sensitivity of the bolometers at 100 mK is limited by 
the photon noise and their short time response ranges from 5 
to 14 ms, which is adequate for the Archeops scanning and 
acquisition strategies. 

2.4. Optical configuration 

The Archeops optical configuration consists of a 1.5 m ofF-axis 
Gregorian telescope illuminating a set of back-to-back horns 
which are coupled to each of the detectors in the focal plane. 
The horns, which are corrugated and flared, are cooled down 
to 10 K by helium vapors and their wave guide sections act as 
frequency high-pass filters. Low-pass filters are located at the 
back of the horns on the 1.6 K stage. These two sets of filters 
define the frequency band of operation of Archeops from 100 
to 600 GHz. The complete set of bolometer, filters and horn 
constitute a photometric pixel. Those are layed out on con- 
stant elevation (scan) lines. The telescope images these lines 
into curved lines on the focal plane. Figure |2] shows the focal 
plane layout of the photometric pixels. The layout of bolome- 
ters at different frequencies was chosen so as to have redundan- 
cies on diff'erent angular and time scales. The main axis of each 
photometric pixel is pointed at the image of the primary mirror 
through the secondary mirror. Entrances of the 10 K back-to- 



Fig. 3. Spectral transmission of the various types of photomet- 
ric pixels. This was obtained by combining different measure- 
ments at the component level. 



Fig. 4. From top to bottom, evolution of the temperature of the 
focal plane and of the 1.6 K and 10 K cryogenic stages during 
the KS3 flight 



back horns are located above the focal plane at various heights 
(typically about 6 mm), in order to prevent optical cross-talks 
between channels. 

2.5. Observation frequency bands 

The Archeops data are acquired at four frequency bands cen- 
tered at 143 (8 photometers), 217 (8 photometers), 353 (6 pho- 
tometers) and 545 (1 photometer) GHz. The first two were ded- 
icated to the measurement of the angular power spectrum of the 
CMB temperature anisotropics. The last two were designed to 
measure the dust diff'use Galactic emission. Figure|3lshows the 
spectral transmission of the various type of photometric pixels 
corresponding to the four Archeops frequency bands. 

The 353 GHz photometers are arranged in pairs cou- 
pled to the same horn via an Ortho Mode Transducer (OMT, 
Chattopadhyay et al. (1999)1 and are optimized to measure the 
polarized sky signal as described in |Benoit et al. (2004), 

2.6. Cryostat monitoring 

Thermometers are used to monitor the cryostat temperature at 
each of the thermal stages described above. The thermometers, 
made of large thin films of NbSi, are added to the 100 mK, 1 .6 
and 10 K stages. They are described in Camus et al. (2000)| 
The house keeping data obtained from these thermometers 
are essential for the subtraction of low frequency drifts in the 
Archeops bolometer data as described in Sect.0 

3. In-flight performance 

We discuss in this section the performance of the Archeops in- 
strument during the KS3 flight in terms of cryogenics and pho- 
tometry. The flight took place on February 7''' 2002 and lasted 
21.5 hours, starting at 12h44 UT time. During 19 hours the bal- 
loon was at nominal altitude, between 32.5 and 34.5 km above 
sea level, with a total of 12.5 hours of night data used for sci- 
entific purposes. 

3.1. Cryogenic performance 

The cryostat functioned autonomously during the entire flight 
duration. The dilution flow was changed twice. First it was de- 
creased at the beginning of the mission to increase the life time 
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of the dilution. Then, it was increased at sunrise to compensate 
for the extra thermal power from the sun. 

The full cryogenic system warmed up mechanically at 
launch, being at nominal temperature at about 15hOO UT time. 
Fig.|3]shows the temperature of the focal plane, the 1.6 K and 
the 10 K stages from top to bottom respectively. At float alti- 
tude, the focal plane cooled down staying below 100 mK dur- 
ing the entire flight. A plateau of about 90 mK was reached 
since 19h00 UT. The 1.6 K stage was stable during the entire 
flight at a temperature of about 1.5 K. The 10 K stage remained 
at about 9 K until sunrise at 27h00 UT^, and then increased up 
to about 12 K. During the night flight the temperature of the 
bolometer bath was stable at 90 mK. 



3.2. Bolometer signal and noise contributions 

3.2.1 . A simple photometric model 

In order to evaluate a priori the total background incoming 
onto each bolometer, we performed a component by compo- 
nent photometric analysis by dealing with emission processes 
and various transmission coefficients. From the sky to the de- 
tector, it includes 

- the CMB emission assuming a simple 2.725 K blackbody 
(Mather et al. (1999) I 

- the atmospheric emission for which the emissivity was 
computed for 41 deg. elevation (airmass of 1.52) at 
32 km altitude using the Pardo's atmospheric model 
([Pardo et al. , (2002)). A temperature of 250 K is assumed. 

- the radiation from the telescope which is assumed to 
have an emissivity of 0.00285 x 2 x VI mm//l 
( |Bock et al. (19 95) l. The factor 2 is to account for the pri- 
mary and secondary mirrors. A temperature of 250 K is 
assumed. 

- the emission of the polypropylene window, which allows 
the radiation to propagate to the cold optics while main- 
taining the vacuum, is neglected here. 

- the radiation from the 10 K stage which was found to emit 
in the KS 1 flight a detectable fraction of the background but 
which was negligible for the KS3 flight. We also account 
for a transmission factor across the lOK stage estimated to 
be 0.6. 

- the transmission curves shown in Fig.|3l They mostly rep- 
resents the filtering done at the 1.6 K stage, although some 
filters were sometimes placed on the 10 K stage and a band- 
pass filter at the entrance of the 100 mK horn. 

- the bolometer with an assumed perfect absorption. 

3.2.2. A simple bolometer model 

The theory of the thermodynamical and electrical behavior 
of bolometers is described in details by .Mather (1984 ) and 
|Piat et al. (200 1)| Here we concentrate on the main equations 
to introduce the parameterization given in subsequent tables. 



We characterize the thermistor behavior of thermometers 
and bolometers using 



R- Roc exp 



(1) 



The electron-photon decoupling is neglected as well as electric 
field effects in the bolometer. The absolute temperature Ti of 
the thermistor is calibrated with carbon resistance and is valid 
within 3 milliKelvin. 

The heat equilibrium for the bolometer reads: 



Pc{Ti,To)^g 



To 



= P, + Pr 



(2) 



where is the base plate temperature measured by the 
standard thermometer, T\ is the bolometer temperature ob- 
tained through R{T\) (Eq. [0, T\qq is a reference temperature 
(we take it as T\qq - 100 mK) so that the constant g is in unit 
of pW, Pc is the cooling power, Pj - UI is the Joule power 
dissipated in the thermistor and Pr is the absorbed part of the 
incident radiative power. 



Bolometer 




r,(K) 


G(pW/K) 




C(pJ/K) 


143K01 


11.27 


21.16 


57.96 


2.80 


0.40 


143K03 


44.23 


16.60 


70.62 


2.30 


0.42 


143K04 


52.91 


16.85 


63.83 


2.30 


0.47 


143K05 


53.41 


16.85 


50.94 


2.30 


0.35 


143K07 


21.55 


18.79 


60.30 


2.55 


0.38 


143X01 


21.47 


20.03 


116.29 


2.85 


0.57 


217K01 


299.42 


16.73 


28.93 


2.10 


0.18 


217K02 


189.54 


13.21 


65.15 


2.00 


0.36 


217K03 


242.93 


13.10 


62.87 


2.00 


0.03 


217K04 


159.03 


13.76 


69.45 


2.20 


0.39 


217K05 


1172.80 


9.72 


59.31 


1.65 


0.34 


217K06 


120.92 


14.18 


65.11 


2.10 


0.46 


217T04 


52.38 


14.79 


161.00 


2.30 


1.06 


217T06 


136.69 


13.67 


182.16 


2.30 


0.00 


353K01 


99.45 


14.39 


90.50 


2.20 


0.21 


353K02 


94.59 


15.01 


98.47 


2.20 


0.12 


353K03 


85.27 


14.86 


99.09 


2.20 


0.22 


353K04 


68.05 


15.18 


106.59 


2.20 


0.23 


353K05 


77.21 


14.67 


103.89 


2.20 


0.39 


353K06 


50.70 


18.49 


116.00 


2.20 


0.38 


545K01 


34.94 


18.48 


136.70 


2.20 


0.04 



Meaning 3hOO UT the next day 



Table 1. Bolometer model parameters as described in Eq.^and 
lllfor each of the Archeops bolometers. 



Table ^ lists the parameters of all thermometers and 
bolometers that were used during KS3 flight (bolometers or- 
dered by channel, 217K05 was blind during the flight). The 
differential conductivity at 100 mK is G = dPc/dT = /Sg/Tim 
(from Eq. These constants are consistent with those mea- 
sured on cosmic rays. A more detailed description of the pre- 
vious issues is given in Sect. |5] The heat capacity was sim- 
ply taken as C = tjG where is the first time constant of 
the bolometer. Time constants are derived from a fit on Jupiter 
data taken during the flight, including electronic filtering and a 
Gaussian beam (see Sect. 16. 1> . 
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The differential conductivity G is taken at lOOmK. Note 
that the "Kiruna" bolometers have typical conductivity be- 
tween 60 and 80 pW/K and heat capacity of 0.3 to 1 pJ/K al- 
though some of them deviate significantly from this range. 

3.3. Detector noise 

We include in the detector noise model contributions from the 
FET electronics, the Johnson noise and the bolometer thermo- 
dynamic noise (.Mather (1984). i. To the detector noise we add 
quadratically the photon noise deduced from our photometric 
model. Note that the bolometer noise is not white at high fre- 
quency due to the bolometer time response. Figure [S] shows 
an example of power spectra of the time ordered data of four 
representative Archeops bolometers during the KS3 flight. We 
overplot the noise model discussed above which is in quali- 
tative agreement without any parameter tuning at frequencies 
higher than a few Hz. We observe an increase of power with 
decreasing frequency mainly due to the low frequency system- 
atics. Although we have smoothed out the power spectrum we 
still can observe peaks which correspond to the sky signal at 
the spin frequency harmonics which are mainly dominated by 
Galactic and atmospheric emissions. Centered at 1 Hz and in 
particular for the high frequency channels there is a very pe- 
culiar structure which may be of atmospheric origin. Finally at 
high frequency we observe correlated structures. A more de- 
tailed description of systematics and their subtraction is given 
in Sect. 121 

Finally, we present in table|2a summary of the noise prop- 
erties of all Archeops bolometers. From left to right we in- 
clude representative values within the night flight for the main 
photometric quantities of those bolometers: current, resistance, 
and responsivity. Next we provide the expected absorbed power 
from a simple photometric model made of subsystem transmis- 
sion measurements presented in section lT.2.11 Then we give the 
absorbed power as measured (with 1 pW absolute uncertainty) 
using the bolometer model described in section 13. 2. 21 The ef- 
ficiency is given as the ratio of Jupiter inflight calibration with 
calibration from the photometric model. Noise measurements 
are given at the bolometer level both for photon and total noise. 
The photon noise is within a factor 2 of the total noise as mea- 
sured in flight conditions. 

4. Preprocessing 

In this section, we describe the preprocessing of the Archeops 
data. This includes the demodulation of the raw data and the 
removal of the parasitic signal introduced by the readout elec- 
tronic noise. We also describe the linearity correction of the 
bolometers and the flagging of the data aflected by cosmic rays 
and noise bursts. 

4.1. Prefiltering 

The data are acquired in total power mode via an AC square 
wave modulated bias. All the modulations are driven by the 
same clock at 76.3 Hz, leading to an acquisition frequency 
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Fig. 6. Top : Kernel of the digital filter used for demodulation 
(see text for details). Bottom : Fourier power spectrum of the 
digital filter compared to a square filter, to the beam pattern and 
to the bolometer time constant. 



face/ - 152.6 Hz. The AC square wave modulated bias trans- 
forms the data into a series of alternative positive and negative 
values. This induces a peak at the Nyquist frequency, facql'^, in 
the Fourier power spectrum of the bolometer data. This peak 
fully dominates the signal and needs to be removed for de- 
modulation. This is performed by filtering the data and for this 
purpose we have constructed a digital filter with the following 
constraints: 

- the transition after the cut-off' frequency, taken to be 60 Hz, 
must be sharp for a complete removal of the modulation 
signal, 

- the ringing of the Fourier representation of the filter above 
the cut-off" frequency needs to be below the 2% level, to 
avoid any possible aliasing. 

These two constraints lead to a digital filter of 23 points 
whose kernel is shown on the top panel of Fig. 14. II The Fourier 
response of the filter is shown on the bottom panel and com- 
pared to that of a simple square filter. We observe that the high 
frequency cutting of the digital filter is much sharper than for 
the square filter therefore preserving the signal better. For com- 
parison, we also plot the Fourier response of the bolometer time 
response and of the beam pattern which determine the spectral 
band for the signal. No signal is therefore removed by the dig- 
ital filter. 

4.2. Removal of readout digital noise 

As discussed above, the bolometer signal is biased with a 
square signal. The data are then amplified by a digital pream- 
plifier and buffered and compressed by the on-board computer 
into blocks. The blocks are then recorded. The compression 
procedure preserves most of the signal of interest. Code 32- 
bit words at the beginning and end of each block allows us to 
check those blocks. These blocks are of different sizes depend- 
ing on the nature of data which may correspond to the signal 
from the gyroscopes, the bolometers, the thermometers or the 
stellar sensor. The length of the bolometric and thermometric 
blocks is of 72 samples. 

Fig. 7. Top: Fourier power spectrum of KS3 143K01 bolome- 
ter data showing the frequency peaks produced by the readout 
electronic noise. Bottom: Same after preprocessing. The am- 
plitude of the peaks is significantly reduced. 



While the on-board computer deals with in-flight com- 
mands, the data recording is delayed and a few data blocks 
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Fig. 5. From top to bottom and from left to right are shown the power spectra (in 10 '^W. Hz 5) of the Archeops 143K03, 
217K06, 353K01 and 545K01 bolometers, respectively, as a function of frequency (in Hz). For comparison we overplot (smooth 
curve) the detector noise contribution for each bolometer as given by the model presented in the text. 
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Table 2. For all Archeops bolometers from left to right. Photometric quantities as representative of night flight values: the current, 
the resistance, and the responsivity. Expected absorbed power from a simple photometric model made of subsystem transmission 
measurements presented in section 13.2.11 Absorbed power as measured (with 1 pW absolute uncertainty) with the bolometer 
model described in section 13.2.21 Efficiency as the ratio of Jupiter inflight calibration with calibration from the photometric 
model. Photon and total noise. 



are buffered before recording. Small offset variations in the 
electronics lead to significant differences between the mean 
value of the last recorded blocks and the next ones. As in-flight 
commands are sent and received by the on-board computer 
periodically during the flight (every few data blocks), the 
differences in the mean between blocks induce a parasitic 
signal on the data. This parasitic signal shows up in the data as 
a periodic pattern of basic frequency facqlT^- Further, as series 
of blocks are buffered before recording we also observe in the 
data periodic patterns at frequencies which are submultiples 
of /flcfy/72. For most of the bolometers this systematic signal 
dominates over the noise and is clearly visible both in the 
time and frequency domain. The top panel of Fig. shows a 
zoom-up of the power spectrum of the data of the KS3 143K01 
Archeops bolometer. We observe on the spectrum a series of 
peaks which correspond to the parasitic signal. 

The subtraction of the parasitic periodic signal can be easily 
achieved using a time domain template for it. Indeed, we have 
implemented a fast algorithm for calculating a time varying 
template of the parasitic signal. First of all, for each Archeops 
timeline we have divided the data into pieces of N blocks of 
720 samples. The block size corresponds to the largest period 
between two in-flight commands. Then each piece of data has 



been reordered into a 720 xA^ matrix so that a time evolving pat- 
tern of the parasitic signal over 720 samples can be calculated 
by smoothing up over the N blocks. The exact number of 720- 
samples blocks to be summed up is a compromise between, 
first, the minimum signal to noise ratio needed for extracting 
the parasitic signal from the data; second, the time evolution 
rate of the parasitic signal and third, the minimum time inter- 
val needed to consider that the sky signal varies sufficiently for 
not contributing to the template. We have found that for most 
Archeops bolometers = 100 is a good compromise. The con- 
structed template is repeated times (size in samples of the 
time interval processed) and then subtracted from each piece 
of data. 



The bottom panel of Fig. shows the power spectrum of 
the KS3 143K01 bolometer after applying the above procedure. 
The procedure reduces significantly the peaks. For example the 
fundamental frequency peak at 2. 12 Hz is reduced to much less 
than 10 % of its original value. The peak at 12.7 Hz, although 
significantly reduced, is still visible in the preprocessed spec- 
trum. It will be cut off in the Fourier domain as discussed in 
Sect.O 
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4.3. Linearity correction 

The cryostat temperature underwent a slow decrease during the 
flight, leading to a slow change of the calibration in mK/fiW. 
This change in calibration can be corrected for by modeling 
the responsivity of the bolometer. Actually, for a TOI b in juV, 
we can write the linearity corrected TOI as follows 



-Vb\n + Vo 

Vh 



(3) 



The parameters Vi, and Vq are determined from the bias- 
current curves of each of the bolometers. After this smooth 
correction, the calibration factor in mK/fiV can be considered 
as constant over the flight, thus allowing for a much easier de- 
termination. For the KS3 flight the correction does not exceed 
20 %, and it is only important for the first 2 hours of flight. 
This has been cross-checked via the Galactic plane calibration 
method described in Sect. 19.21 

4.4. Flagging of the data 

In this section, we describe the identification and flagging of 
parasitic effects including glitches, noise bursts and jumps in 
the data. For Archeops most glitches are due to the increase 
of temperature of the bolometer due to the energy deposited 
by cosmic rays hits. Jumps are essentially due to changes of 
the equilibrium voltage of the bolometer and there are only 
a few during the whole flight. Bursts of noise are induced 
by microphonic noise coming mainly from the mechanical 
oscillations of the gondola. 

To flag and remove the data affected by the above system- 
atic effects, the first step is to detect spikes in the TOIs above 
a certain threshold level. For this purpose the r.m.s. noise level, 
cr^, is estimated locally on a 400 points running window as the 
standard deviation from the median value, m, of the data after 
removing 5% of the lowest and largest data values. The data 
with flux above 8cr are considered as glitches. To preserve the 
Galactic signal, which can be sometimes spiky or/and larger 
than the threshold limit, a baseline fbase fitted as a combina- 
tion of the two first Fourier modes is removed whenever data 
values above the threshold are detected at Galactic latitude be- 
tween -10° and H-10°. The value of cr is then re-computed and 
the above criteria re-applied. This technique is time-consuming 
but not required outside of the Galactic plane where a flat base- 
line is already a very good approximation. 

The second step is then the flagging of the data. When the 
parasitic signal is due to a glitch, we can model it by the con- 
volution of a Dirac delta function at time f, with the sampling 
window and a double decreasing exponential function with two 
time constants t short and Tio„g. The first corresponds to the re- 
laxation time of the bolometer itself and the second is of un- 
known origin. The time constant values depend only on the 
bolometer and must be the same for all glitches impacting this 
bolometer The main objective of this first analysis is not to re- 
produce faithfully the glitch shape but to estimate which part of 
the data is badly affected by it and must be flagged. Therefore, 
the same conservative values are adopted for all bolometers. 



T short - 2 samples (13 ms) and Tio„g = 50 samples (325 ms). 
We then fit the following glitch model 



/(f,f,) 



^ short * 



facq + ft 



base 



(4) 



where the free parameters are the amplitudes of exponential 
functions A s/,o,, and A;o„j and the baseline and where facq is the 
sampling frequency. All data samples within f„„„ and fmoi given 
by 



= f ,■ - 11 

= h + T short In 



O.lo- 



+ T long In 



^long 



O.lcr 



-1-11 



are flagged. So data samples for which the glitch contri- 
bution is at a level higher than 10 % of the local noise are 
flagged. The extra 1 1 samples margin on each side of the glitch 
position accounts for the effect of the digital filter at 23 points 
discussed in the previous subsection. An additional margin of 
100 samples is used when the fits is of poor quality. This extra 
flagging concerns about 20 % of the glitches detected on the 
OMT bolometers, 33 % on the Trapani-like bolometers and 
less than 15 % for the others bolometers of the KS3 flight. 

Detailed statistics of the number of glitches detected in the 
KS3 flight are reported in Tab. |3l The bolometers from the 
Trapani flight are quite sensitives to glitches, 15 to 20 glitches 
per min. Polarized OMT at 353 GHz shows a rate of ~4 glitches 
per minute, whereas at 217 GHz and 545 GHz we detect less 
than 2 glitches per minute. At 143 GHz the bolometers present 
a glitch rate between 1.5 and 4 per minute. The glitch rate is 
fully related to the effective surface of the bolometer which 
varies between bolometers. A larger glitch rate can be ex- 
plained by a larger effective area of the spider-web absorber. 

In the above procedure bursts of noise on the data are as- 
similated to glitches and the flagging obtained is poor To en- 
sure a better flagging, we proceed to a visual inspection of the 
data. We check all the pieces of data found above the thresh- 
old limit and extend manually the flagging if necessary. Those 
data samples affected by noise bursts are flagged as such. We 
also observe jumps on the data which are caused by extreme 
changes on the DAC currents of the bolometer. 

The values of the DAC currents are stored as housekeeping 
data and allows us to correct the data from those jumps via a 
simple destriping algorithm. In addition, by visual inspection 
we determine the data samples which are affected by jumps 
and they are manually flagged. 

At the end of the process, a total of 1-2.5 % (resp. 2-4 % 
and 12-18 %) of the data are flagged for the KS3-like bolome- 
ters (resp. for the OMT and Trapani-like bolometers). Flagged 
data are then replaced by constrained realization of noise as 
discussed in Subsect.l84l 



5. Pointing reconstruction 

The knowledge of the pointing attitude was not needed during 
the flight but its accurate a posteriori reconstruction is critical 
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Bolometer 


# glitches [per min.] 


data flagged [%] 


143K01 


1.8 


0.93 


143K03 


3.6 


1.58 


143K04 


4.2 


1.74 


143K07 


1.6 


0.77 


143K05 


2.2 


0.94 


143T01 


16.8 


8.58 


217K01 


1.0 


0.44 


217K02 


1.1 


0.54 


217K03 


1.3 


0.55 


217K04 


1.6 


0.79 


217K05 


1.3 


0.58 


217K06 


1.5 


0.79 


217T04 


16.9 


8.43 


217T06 


20.7 


11.62 


353K01 


4.8 


2.15 


353K02 


4.1 


1.88 


353K03 


5.7 


2.55 


353K04 


3.8 


1.72 


353K05 


4.9 


2.22 


353K06 


3.3 


1.52 


545K01 


1.1 


0.76 



Table 3. Statistics of glitches per minute for the KS3 flight and 
proportion of flagged data. 



for mapping correctly the sky signal. The pointing of each of 
the detectors in the focal plane is computed as follows. First, 
a pointing solution for the payload is obtained from the pro- 
cessing of optical data collected by the fast stellar sensor (FSS) 
during the flight. Finally, we estimate the pointing offset with 
respect to the payload axis for each bolometer using the recon- 
struction of the focal plane from measurements of point sources 
(see Sect.|6lfor details). 

We have developed an algorithm to extract star candidates 
from the FSS time-sampled photodiode signals (see Sect. 12. 2> . 
Each star candidate is kept into a table including its detection 
time, its position along the diode array and the electrical in- 
tensity observed. The position of the star candidate along the 
diode array is given in terms of an effective diode number. The 
electrical intensity measured is proportional to the logarithmic 
value of the flux of the star. 

5.1. Coordinate system definition 

In the following, we use equatorial coordinates {a, 6} to define 
the position of celestial objects on the sky. The FSS data are 
also easily handled with local coordinates associated with the 
gondola frame, for which the zenith corresponds to the gon- 
dola spin axis direction. The direction of a star on the celestial 
sphere is then given by 6, the angular distance between the spin 
axis and the direction of the star (hereafter the axial distance), 
and by tp, the phase corresponding to the azimuth measured 
from the North. 

To reconstruct the pointing direction of the gondola we 
need to find the direction of the center of the diode array. The 
instantaneous pointing solution is fully described by the set 
{ap, 6p, (f), where Up and dp are the equatorial coordinates of 
the gondola spin axis and (f, the phase for the diode array. Note 



that the phase value ip is the same for all the diodes in the ar- 
ray, therefore also the same for all detected stars, because the 
diode array is placed perpendicularly to the scanning direction. 
In other words, the number of the diode lightens is only given 
by the axial distance of the observed star. 

5.2. Reconstruction Metliod 

The goal is to produce an optimal pointing solution as a func- 
tion of time. The reconstruction is based upon the comparison 
between FSS data and a dedicated star catalog compiled from 
the Hipparcos catalog. The electrical intensity of stars in the 
catalog is computed by taking into account the FSS spectral 
response. Hereafter, we call signal a star candidate in the list 
produced by the FSS software and star an object taken from 
the star catalog. First of all, we find in the star catalog the best 
star to be associated which each FSS signal. We call identi- 
fied signal a signal for which this association is performed. Via 
these associations we obtain a pointing solution for each iden- 
tified signal. Finally, we fit the overall set of identified signals 
through the scan path to get a pointing solution as a function of 
time. 

5.3. Initial pointing estimate 

To be able to associate signals to stars, a first estimate of the 
pointing solution is needed. This is obtained via the GPS data 
which give the local vertical direction, which corresponds to 
the spin axis direction {ap, dp} to an accuracy of ~ 1 degree, 
taking into account the gondola average pendulation. Then, we 
match signal and star directions and try to identify for each 
signal the corresponding star. There is no direct measurement 
of the FSS phase (p. We need to reconstruct it from the rotation 
period by integration. 

5.3.1 . Rotation period 

We now describe the gondola motion relative to the celestial 
sphere. We thus call rotation period the elapsed time between 
two successive detections of the same star after one revolution. 
Each revolution takes about 30 seconds. Due to the Earth mo- 
tion, the spin axis moves about 5' in a per revolution. Each 
star can thus be seen several times by the FSS. For each sig- 
nal, we look for every compatible signal seen in the preceding 
revolution. A compatible signal has a similar intensity and a 
nearby diode number. Time differences between the signal and 
those seen in the last revolution are binned into an histogram. 
The most populated bin gives us the rotation period. Figure |8] 
shows the evolution of the rotation period as a function of time 
for 3 hours of the KS3 flight. This evolution is mainly due to 
the presence of strong stratospheric winds during the flight. 



Fig. 8. Rotation period evolution during the KS 3 flight. 
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5.3.2. Star Sensor Phase 

We reconstruct the FSS phase by integrating the angular speed 
where T is the rotation period. The resulting estimate (p dif- 
fers from the phase by a slowly varying offset. To correct from 
this bias, we compare, for each revolution, the phase of the 
most intense signals with the phase of the brightest stars lo- 
cated in the 1 .4 degree wide band scanned by the diode array 
during a revolution. 

Fig. 9. Evolution of the distribution of phase differences be- 
tween signals and bright stars for the KS 3 flight. 



Fig. 10. Distribution of the axial distance of bright stars versus 
the diode number of the corresponding intense signals. Notice 
the strange behavior of the diode 26. This diode is excluded 
from analysis. 



The analysis of phase differences if*-(p gives us the FSS 
phase offset shown in Fig.|9l The distribution of axial distance, 
9, values of bright stars associated to intense signals allows us 
to adjust the geometrical relation between the axial distance 
and the diode coordinates along the array. Figure^] shows the 
distribution of the axial distance of bright stars as a function 
of the diode number of the corresponding intense signal in the 
FSS. We observe for each diode number the distribution has 
a well defined peak from which we can reconstruct the axial 
distance for each diode. The width of the peak is due to the 
pendulation motion of the gondola. 

5.4. Star-signal matching algorithm 

The association algorithm used above is based on a comparison 
of the star and signal directions. An error 6 on the spin axis 
direction (ap, 6p) translates into a local rotation and then an 
error on the reconstructed direction for each of the signals. The 
gondola pendulation is a slow time varying function on scales 
of a few degrees. Therefore, the error d and the local rotation 
parameters, are slowly varying functions too. In other words, 
signals detected within a few degrees area are thus shifted a 
roughly equal amount from their true position on the celestial 
sphere. 

The matching algorithm is based on the above statement 
and proceed as follows. First, for each signal, we associate stars 
and signals with compatible positions and intensities. Second, 
given a reference signal, we check whether for the following 
signals there are N stars such that the corresponding shifts are 
close. If so, this displacement is the signature of a local rotation 
induced by a wrong reconstruction of the spin axis direction or 
by a wrong estimate of the FSS phase. 

Free parameters like the number N of signals used or the 
tolerance on the angular distance between the signal and its 



corresponding star directions can be tuned to optimize the as- 
sociation efficiency. In practice, tight cuts on those parameters 
reduce the probability of wrong associations but at the same 
time reduce the number of good associations available on the 
regions where the pointing reconstruction is bad. To improve 
this situation we use the fact that the FSS sees a given star dur- 
ing several revolutions. Once a good association is obtained 
we propagate this information to the whole data set using our 
estimate of the rotation period and thus we can improve the 
association efficiency and therefore the pointing solution. 

5.5. Pointing solution improvement 

The axial distance 6 is the only quantity which can be directly 
measured. When signals have been associated to catalog stars, 
it gives a way to reconstruct the spin axis direction. As the po- 
sition of the signal and that of its associated star must be the 
same, the spin axis is therefore located on a cone centered on 
the star with an opening angle 9. Using two couples, signal- 
star, we can find the direction of the spin axis. Indeed the in- 
tersection of the two cones, one for each couple, leads to two 
solutions. Only one of them is geometrically relevant. Using 
the whole data set we can thus correct the estimate of the spin 
axis direction during the flight. We upgrade the FSS phase tak- 
ing into account the new estimate of the spin axis pointing. The 
process is iterative to obtain a more accurate estimate of the 
pointing for the whole flight. The increase in accuracy at each 
iteration allows us to use tighter cuts to get a better quality 
matching between stars and signals. 

The FSS dataset available is mainly composed of faint stars 
making the above iterative solution very important. Further, 
calibration uncertainties on the signal get broader as the inten- 
sity decreases. The associations for the brightest stars allows us 
to recalibrate the FSS signals. Adding finer constraints on the 
intensity of the signal increases the quality of its association 
to a star in the catalog. This also improves the final pointing 
solution. 

5.6. Scan path fit 

Once the signal-star associations are obtained we have a dis- 
crete pointing solution at the times where the signal were ob- 
served. Our purpose is to generate an optimal continuous scan 
path and then we have to interpolate the pointing solution along 
FSS data. This solution should not only be interpolated, but 
also optimized along the data set {jc,, f,). To get an optimal 
pointing solution we have to reconstruct {ap(t), 6p(t), ip(t)), 
from the {apt, dpi, (pi}, set, where / labels a given signal-star 
association. This is performed by computing first a smooth so- 
lution for the pointing and then correcting it. 

Smooth pointing solution 

We first produce a smooth solution for the pointing (a°, 6^^ and 
by fitting the set {a^/, 6pi, tfi}, using a chi-square mini- 
mization. As the set [apt, dpi, (pi}i is irregularly sampled in 
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time we obtain a generic interpolation, x(t), of the pointing so- 
lution through a decomposition of the form 



x(t) ^Y^CkU{t-tk) 



(5) 



Fig. 11. 2D histogram showin'g a variation of phases differ- 
ence between signals and associated stars with the signal diode 
number. 



where each U{t - tk) is a generic kernel of the form 



this direction induces a modification of the coordinates 9* and 



sinc{^)e -p-"^ centered alt - tk- We choose the parameter pio Xo first order, we have 
optimize the representation of the low frequency components 
in x{t) and cr. 

The coefficients {q) are obtained from the minimization of 
the chi- square 

leading to the following linear system 

N ( \ 

^ Y^u(ti-ii}U{ti-tk} ck 

k=\ V i ) 

= J^XiU(ti-ii) with k,l^ 1,2,... ,N. 

i 

We solve this system for the three quantities of interest 
ap(t), 6p(t), and ip{t). 



if* = ip*° + C2i^ap + C22Mp. 



The coefficients c,y are known functions of ap°, 6p", 0*° and 



Then Eq.|6lbecomes 



2 _ \ ^ {Vf -V>° + C2i^ap + C22A6p - Af 



X 



9* ° - 6i + ciiA&p + cuASp 



Corrected pointing solution 

Once we have a first smooth solution for the pointing (5°, 5° 
and ifP) we compute corrections to it Aap(t), Adpit) and Aip(t). 
For this purpose we decompose these 3 quantities in terms of 
kernel functions as in Eq.|5l We call Aait, Adk and Apk the de- 
composition coefficients for A&pif), A6p(f) and Aip{t) respec- 
tively. In this case we consider high frequency terms to opti- 
mize the pointing solution. 

The FSS dataset [Upi, dpi, (pi}i can be rewritten more ex- 
plicitly as {up i, 6p i, (fi, 6i, a*, 5*),-. 0, is a linear function of the 
diode number a* and 6* are the coordinates of the star corre- 
sponding to signal /. This set can also be expressed for the star 
position in gondola frame coordinates {up ,, dp (pi, 9i, (f*, 0*),-. 
We can obtain an estimate of the pointing corrections by com- 
paring the reconstructed star positions with the pointing posi- 
tion at the time of their observation. 



To achieve convergence, we follow the same steps as the 
method described in l5.5l We first calculate the correction on the 
spin axis direction Aap(t) and Adp(t). This leads us to minimize 
the quantity 



2 _M-^('P*" -ip° + C21 Affp + C22A6p 



3*° -Oi + cuAap + cnASp 



(7) 
(8) 



The phase correction is then obtained by taking 



2 sr\^i "^zi A&p + C22A6p - Aip 

X2- 



~\2 



(9) 



X 




(6) 



We note cr* and o"*' the errors associated with the measure- 
ments of and (fi obtained in the previous section. There are 
two sources of asymmetry between cr^' and cr^. The axial dis- 
tance coordinate is directly measured by the FSS. The phase 
coordinate must be reconstructed, once the spin axis direction 
is known. The second source is geometric. A diode covers 1.9 
arcminutes in the cross-scan direction by 7.6 arcminutes along 
the scan. We use cr^ = cr and erf = 2cr. 

The star coordinates (p* and 0* in the gondola frame de- 
pend on the spin axis direction. A variation A&p and A6p in 



The minimization of X2 ^^'^ X2 l^ads to the iterative 
resolution of two linear systems with free parameters Aat, Adk 
and Apk- 

In the above we have assumed that the photodiodes array 
was oriented perpendicularly to the pointing direction. This hy- 
pothesis can be verified by comparing the phase for the stars 
with the phase of the FSS as a function of the diode number. 
This comparison is shown on Fig.[^ We observe a phase shift 
which indicates that the photodiodes array is tilted along the 
scan direction. Given the 1.8 m focal length of the parabolic 
mirror and a 1mm photodiode area along cross-can direction, 
we find an inclination of ~ 3 degrees. The phase of each signal 
is thus corrected to take this effect into account. 
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5.7. Accuracy 

We have of two independent but complementary ways of as- 
sessing the accuracy of the Archeops pointing reconstruction. 
A first estimate can be obtained from the distribution of co- 
ordinate differences between the signals and their associated 
stars. Figure[2]show the distribution of errors in the plane ax- 
ial distance-phase before and after scan path fit respectively. 
The 95% and 68% confidence level contours are displayed in 
white. We observe that the axial distance coordinate has intrin- 
sically a better accuracy by a factor 2.5. Further, we notice a 
significant decrease of the errors for both the axial distance and 
the phase. 

We can also estimate the errors in the pointing reconstruc- 
tion via the Fisher matrix of the free parameters in the scan- 
path fit described by Eq|Hl and |5] This gives us a continuous 
estimate of the pointing error which is used to flag those re- 
gions for which the pointing is badly known. Hereafter, we call 
this flag on the data bad pointing flag. The distribution of equa- 
torial coordinate diff'erences in Fig.[0]shows that the attitude 
reconstruction is achieved with an accuracy better than 1 .5 and 
1 arcmin in RA and DEC respectively, at the 1-cr level. 

Fig. 12. From top to bottom, distribution of errors in axial dis- 
tance - phase plane with 95% and 68% confidence levels (in 
white) before and after scan path fit respectively. 



observations obtained at 15h36 and 18h427 UT hours to cross 
check the results. Saturn can be also considered as a point 
source at the Archeops resolution (apparent diameter of 19 by 
17 arcsec). 

For each detector, we start by computing local maps of the 
planets in azimuth-elevation coordinates which correspond to 
the along-scan and cross-scan directions. These maps are ob- 
tained by projecting the TOl data without filtering. To remove 
the zero level in these maps we estimate a baseline in the TOl 
which is then subtracted. The latter is estimated from a TOl 
where all the flagged data are interpolated using a constrained 
realization of noise. The TOl signal obtained for planet obser- 
vations is the superposition of two main effects. First, the con- 
volution of the source sky signal with the beam pattern of the 
photometric pixels. Second, the convolution of the bolometric 
signal with the time response of the bolometers which is char- 
acterized by a time constant. Both eff'ect are clearly visible in 
the Jupiter map shown on the top panel of Fig. The beam 
pattern convolution widens up the point source signal both in 
the cross-scan and along-scan directions. The effect of the time 
response convolution appears as a tail in the map along the scan 
direction. 

In our analysis we first estimate the bolometer time con- 
stants for each using in-scan profiles of the Jupiter or Saturn 
map. Then we deconvolve the TOl from the bolometer time 
constant and recompute local maps as the one presented on the 
bottom panel of Fig. [H] From these maps, we characterize the 
beam pattern of the photometric pixels. 



Fig. 13. KS3 flight 95% end 68% confidence levels for eiTor 
distribution in equatorial coordinates after scan path fit. 



6. Bolometer response 



Fig. 14. Top: map of Jupiter for the 143K03 bolometer in yuV 
before time constant deconvolution. The map is represented in 
the par and cross scan direction in arcmin. Bottom: As above 
but after deconvolution from the time constant. 



6.1. Time response 

6.1 .1 . Optical time constants estimate 

The time response, TR, of the bolometers can be described by 
the combination of two decreasing exponentials with time con- 
stants Ti and T2 



TR(t) ^(l-a)e 



-t/Tl 



(10) 



with a a mixing coefficient to be estimated from the data. 
As Archeops scans the sky at roughly constant speed, the 
eff'ect of beam pattern and time response are degenerate in the 
along-scan direction. In order to have the simplest possible 
model, but that allows us to separate both effects, we will 
assume that the beam pattern shape is symmetric along the 
scanning direction. ^ 



We describe in this section the reconstruction of the 
Archeops focal plane parameters for the KS3 flight. For this 
purpose we estimate the time response of the bolometers, 
the optical response of the photometric pixels and the focal 
plane geometry on the celestial sphere. The focal plane is 
reconstructed using planets observations. The brightest one, 
Jupiter, is observed twice at ~ 16.5 and 21h00 UT hours and 
can be considered as a point source at the Archeops resolution 
(apparent diameter of 45 by 42 arcsec). We also use Saturn 



The time constants are estimated fitting the Jupiter pro- 
files using a minimization for a grid of 3 parameters ti, 
T2 and a which are set in the range [1,10] ms, [10,100] ms 
and [0,1] respectively. The profiles used are the 4 arcmin cross- 
scan average of local maps of the two Jupiter crossings. For 
each set of parameters, we deconvolve the initial TOl from the 
TR(Ti,T2,a). We then fit a Gaussian on the rising part of the 
beam profile and use a Gaussian with the same amplitude and 

^ Notice that for a constant angular speed the effect of the time re- 
sponse convolution can be assimilated to a beam convolution. 
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sigma for the decreasing part. We do not account for the higher 
part of the beam shape that can present several maxima (es- 
pecially multimode ones) by fitting only the lower 80% of the 
profile data. 

We compute the minimum of the;\f^ in the (ti, T2), (ri, a) 
and {t2, a) planes. The best-fit parameter values are obtained 
from the average of the two maxima obtained for each param- 
eter By integrating the surface we obtain directly the Icr error 
bars. If a is compatible with within Icr, we re-compute the 
estimation of ti using a single time constant model to reduce 
the error bars. 

Figure [21 shows one of the Jupiter map profile for the 
217K04 bolometer before (in red) and after (in black) deconvo- 
lution from the bolometer time response. The tail in the profile 
due to the time response is clearly suppressed after deconvolu- 
tion. 



Fig. 15. 217K04 Beam profile on Jupiter before (in red) and 
after (in black) deconvolution of the two time constants (tj = 
5.5T_\ °l ms, a = OAS^^'-^ and = 39,3S+J/^^ ms). 
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Table 4. Bolometer time constants for the Archeops KS3 flight 



In Tab. 13 we present, for each of the Archeops bolometers, 
the values of tj (ms), a and T2(ms) obtained from the analysis of 
the Jupiter profiles. The analysis of the Saturn profiles provides 
consistent results. 



6.1 .2. Estimate of the bolometer time constant from 
glitclies 

The time response of the bolometer can also be estimated us- 
ing the signal from cosmic ray glitches with short time constant 
(see Sect. 14. 4> which hit directly the bolometer. To a very good 
approximation, for these glitches, the signal is just the convolu- 
tion of a Dirac delta function with the bolometer time response 
and the sampling kernel and therefore they have all the same 
shape. A template of this can be obtained by piling up all short 
glitches in the data after common renormalization both in po- 
sition and amplitude. We can then fit to this template the glitch 
model (Eq.0} to estimate Tshon- Notice that only a few glitches 
have an additional significant long time constant preventing us 
from reconstructing an accurate template. 

Figure^Jshows in black the glitch template computed for 
the bolometer 21 7K01. In red we trace the best glitch model fit 
for it, corresponding to a time constant of 7.8 + 0.1 1 ms. 



Fig. 16. Glitch template for the bolometer 217K01. A single 
time constant model have been fitted to the data. For the best 
fit, traced in red, the time constant is 7.8 ± 0.11 ms. 



In Fig. El the values of Ts/,o,, obtained from the fit of 
glitches are compared to the shorter optical time constant ti 
estimated using Jupiter (Tab|4}- These are compatible within 
Icr for a large number of detectors. The observed discrepancies 
may be due to the intrinsic degeneracy between time response 
and beam pattern and the way we break it. Equally, we can 
imagine different time delays in the detector response depend- 
ing on where exactly the glitch hits. 

The second time constants Tiong and T2 differ by at least 
one order of magnitude. The long time constant measured on 
glitches can be interpreted as a longer thermal link for glitches 
which hit the immediate surrounding of the bolometer and will 
not be considered in the following. By contrast, the second op- 
tical time constant, T2, must be taken into account to accurately 
deconvolve from the bolometer time response. 

6.2. Optical response 
6.2.1. Beam pattern model 

After deconvolving the Archeops TOI from the bolometer time 
response, we reconstruct local maps of Jupiter to estimate the 
beam pattern shape. The beam patterns for the Archeops pho- 
tometric pixels happen to be asymmetric in particular for the 
multimode systems (all the 217 GHz detectors but 217K01, 
217K02 and 217K05 and the two 545GHz detectors). 

We model the main beam shape for each photometric pixel 
using the Asymfast method, described in ,Tristram et al. (2004)| 
This method is based on the decomposition of the main beam 
shape into a linear combination of circular 2D Gaussian. This 
allows us to accurately and simply represent asymmetric beams 
and to convolve full sky maps with them in a reasonable 



J. F. Macias-Perez et al.: Archeops In-flight Performance, Data Processing and Map Making 
Fig. 17. Comparison between glitch and Jupiter short time constant estimates. 



13 



amount of time. This is of great interest when producing sim- 
ulations to estimate the beam transfer function in the spherical 
harmonic plane ( Tristram et al. (2005b)t . 

The Archeops main beams have been modeled using up to 
15 Gaussians. The residuals after subtraction of the model from 
the Jupiter local maps are less than 5%. Figure [Tsl shows an 
example of this multi-Gaussian beam modeling for the photo- 
metric pixel 143K03 using 7 circular 2D Gaussians. From top 
to bottom and from left to right we represent the beam pattern 
shape from the Jupiter local maps, its multi-Gaussian fit, the 
residuals and the histogram of the residuals (black line). The 
latter are shown to be compatible with Gaussian distributed 
noise (red line in the figure). 



Fig. 18. From top to bottom and from left to right, for the photo- 
metric pixel 143K03, the beam pattern map in f/V from Jupiter 
observations, its multi-gaussian decomposition using seven 2D 
Gaussians, the residuals and their distribution which is compat- 
ible with Gaussian distributed noise (red line). 



The resolution for each of the photometric pixels has been 
estimated from a 2D elliptical Gaussian fit to the local beam 
maps from Jupiter observations. The FWHM values given in 
Tableware the geometric mean value in the two directions. We 
present as well the ellipticity as computed from the ratio be- 
tween the minor and the major axis. For monomode systems 
the beamwidth is about 1 1 arcmin. The multimode systems at 
217 GHz and 545 GHz have larger beams with FWHM -15 ar- 
cmin. The 353 GHz beams are monomode but are by construc- 
tion illuminating a small part of the primary mirror to have 
clean polarized beams. A degradation of the beams has been 
noticed in the KS3 flight compared to ground-based measure- 
ments (which are close to the diffraction limit) due to a slight 
out-of-focus position of the secondary mirror after the crash 
landing of the KS2 flight. 

No specific treatment have been applied in the data analysis 
to account for far side lobes as from Planck optical modeUza- 
tion they are expected to be at the percent level. 

6.2.2. Focal plane reconstruction 

The position of each photometric pixel in the focal plane with 
respect to the Focal Plane Center (FPC) is computed using 
Jupiter observations. This then allows us to build the pointing 
of each pixel using the pointing reconstruction as described in 
Sect. Is] In practice we use the relative positions of the 2D cir- 
cular Gaussians of the Asymfast decomposition to estimate the 
center of the beam in focal plane coordinates. Fig.^]shows the 
reconstruction of the Archeops focal plane in azimuth (along 
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Table 5. Resolution in terms of the FWHM and ellipticity of the 
beam pattern for the Archeops photometric pixels in the KS3 
flight. See text for details. The optical mode of the bolometers 
is indicated in the last column: S for single mode bolometers 
and M for multimode bolometers. 



scan) and elevation (across scan) coordinates. The Archeops 
focal plane is about 2 degrees high and 2.5 degrees wide. 

Fig. 19. Focal plane of Archeops reconstructed using Jupiter 
observations. 



7. Description and subtraction of systematics 

In this section we describe in details the main systematic effects 
on the Archeops data as well as the methods and algorithms 
used for their subtraction. Because of the circular Archeops 
scanning strategy the sky signal shows up on the data at fre- 
quencies which are multiples of the spin frequency. This nat- 
urally leads to three distinct regimes in frequency. First, the 
very low frequency components at frequencies well below the 
spin frequency (/ < 0.01 Hz) which are mainly due to 1 //-like 
noise and systematics. Second, the spin frequency components 
(0.03 < / < 3 Hz) at frequencies which contain most of the sky 
signal of interest. And finally the high frequency components 
at frequencies much larger than the spin frequency (/ > 10 Hz) 
which are dominated by detector noise. 
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Fig. 20. Left column: from top to bottom, raw Archeops data in //V (black curve) and reconstructed very low frequency drift 
(red curve) for the 143K03, 217K06, 353K01 and 545K01 bolometers respectively. Right column: from top to bottom, Archeops 
data in juV after removal of the very low frequency parasitic signals for the 143K03, 217K06, 353K01 and 545K01 bolometers 
respectively. 



7.1. Very low frequency systematics 

At very low frequency the Archeops data are dominated by sys- 
tematics coming mainly from temperature fluctuations of the 
three cryogenic stages at 100 mK, 1.6 K and 10 K and from 
the variation of air mass during the flight due to changes in the 
balloon altitude. The left column of Fig. from top to bot- 
tom, shows the raw Archeops data for the 143K03, 217K06, 
353K01 and 545K01 bolometers in the period from 14h00 to 
29h00 UT time. We observe a very low frequency drift in the 
data which is very well correlated within bolometers and also 
with the low frequency components in the housekeeping data 
from thermometers placed at each of the cryogenic stages and 
with measurements of the altitude of the balloon. 

This drift is removed via a decorrelation analysis using as 
templates the housekeeping data described above and a fifth 
order polynomial defined in the time interval of interest. To 
compute the correlation coeflicients we first smear out and un- 
dersample both the Archeops data and the templates, and then, 
we perform a linear regression. A final estimate of the drift 
is obtained from the best-fit linear combination of the tem- 
plates which are previously smoothed down to keep only the 
very low frequency signal. In the left column of Fig. |20| we 
overplot in red the reconstructed very low frequency drift for 
the four Archeops bolometer. In the right column of the fig- 
ure we show the Archeops data after subtraction of the low 
drift estimate which reduces the signal amplitude by a factor 
of 50. Although the decorrelation procedure is very efficient, 
we can still observe a correlated low frequency parasitic signal 
between bolometers. To avoid the mixing up of the bolome- 
ter signals at this stage of the processing of the data, this ef- 
fect is considered within the map making pipeline described in 
Sect.d 

7.2. Spin frequency systematics 

The Archeops data at the few first multiples of the spin fre- 
quency are particularly important. They contain the large an- 
gular scale signal of the sky whose mapping is one of the main 
purposes of the Archeops experiment. At these frequencies in 
addition to the CMB anisotropies and the Galaxy, two main 
components can be identified, the CMB dipole and the atmo- 
spheric emission. Although the former is critical for the calibra- 
tion of the Archeops data at 143 and 217 GHz (see Sect.lgTl 
for this analysis we consider both of them as systematics. 

As above, we perform a decorrelation analysis to re- 
move these parasitic signals. A template for the dipole in 
niKcMS is simulated using the CORE best dipole solution 
( |Smoot et al. (1991)j and the local velocity vector during the 
flight. The atmospheric contamination is mainly due to the vari- 



ation of the air mass induced by changes in the altitude and 
the pointing elevation of the payload. As templates for this ef- 
fect we use the housekeeping data corresponding to measure- 
ments of the altitude of the balloon and the reconstructed el- 
evation of the focal plane. To compute the correlation coeflS- 
cients, we smear out and undersample both the Archeops data 
and the templates and perform a linear regression. The final 
estimate of the spin frequency systematics are obtained from 
the best-fit linear combination of templates. The left column 
of Fig.l^shows from top to bottom, the Archeops data in fjV 
(black curve) and the reconstructed spin frequency systemat- 
ics (red curve) for the 143K03, 217K06, 353K01 and 545K05 
bolometers respectively. We observe that at the low frequency 
channels the dipole contribution (exactly at the spin frequency) 
dominates while at the high frequency the atmospheric emis- 
sion is much more important and dipole becomes negligible. 

The spin frequency systematic estimate is then filtered out 
in the range of frequency of interest, low passed to remove very 
low frequency drifts in the templates and high pass to reduce 
the noise in the templates, and subtracted from the Archeops 
data. The right column of Fig.|^shows the Archeops data for 
the four bolometers before and after subtraction of the spin fre- 
quency systematic estimates. After subtraction we clearly ob- 
serve at all frequency the Galactic signal. This shows up as 
peaks at each revolution and of increasing amplitude with in- 
creasing frequency channel. This is better shown on Fig. |22 
where we represent the power spectrum in yuV yfHz of the 
Archeops time ordered data before (black curve) and after 
(red curve) subtraction of the spin frequency systematics. At 
143 and 217 GHz the dipole contribution appears clearly as a 
peak at the spin frequency before subtraction. After subtraction 
we clearly distinguish the Galactic signal which dominate the 
spectrum at 353 and 545 GHz. 

Beside the atmospheric signal which is correlated to tem- 
plates there is a residual parasitic atmospheric signal. The lat- 
ter can be qualitatively reproduced by simulating turbulent at- 
mospheric layers drifting across the Archeops scanning beam. 
Typical gradients of half a mK^y over 10 degrees azimuth are 
observed at 545 GHz with an evolving period of 1000 seconds. 

7.3. l-ligh frequency systematics 

The Archeops data present, at frequencies / > 5 Hz, parasitic 
noise which shows up on the time domain as noise bursts which 
are neither stationary nor Gaussian. In the Fourier plane this 
parasitic signal is observed in the form of well defined cor- 
related structures or peaks. This parasitic noise is very likely 
related to microphonic noise. Its contribution to the data was 
significantly reduced in the KS3 flight with respect to the KSl 
flight by increasing the distance between the payload and the 



J. F. Macias-Perez et al.: Archeops In-flight Performance, Data Processing and Map Making 



15 



Fig. 21. Left column: from top to bottom, Archeops data in //V (black curve) and reconstructed spin frequency systematics 
(red curve, including the CMB dipole and atmospheric emission) for the 143K03, 217K06, 353K01 and 545K01 bolometers 
respectively. Right column: from top to bottom, Archeops data in /iV before (black curve) and after (red curve) removal of the 
spin frequency systematics for the 143K03, 217K06, 353K01 and 545K01 bolometers respectively. 



Fig. 22. From left to right and from top to bottom power spectrum in /iV ^JUz of the time ordered Archeops data before (black 
curve) and after (red curve) subtraction of the spin frequency systematics for the 143K03, 217K06, 353K01 and 545K01 bolome- 
ters. 



Fig. 23. From left to right and from top to bottom power spectrum in jiW yf{Hz) of the low frequency processed Archeops data 
before (in black) and after (in red) decorrelation from the high frequency noise for the 143K03, 217K06, 353K01 and 545K01 
bolometers respectively. 



spinning motor. The latter was placed far away in the flight 
chain with no rigid link between them. Fig. 1231 shows in black 
the power spectra of the Archeops low frequency processed 
data for the 143K03, 217K06, 353K01 and 545K01 bolome- 
ters. Structures around 15, 22, 40 and 50 Hz are clearly visible 
in the power spectrum. We can also observe a peak at 28.5 Hz 
in the power spectrum of the 353K03 bolometer data. 

These structures are also observed on the power spectrum 
of the data from the focal plane thermometers and from the 
blind bolometer. Moreover, they are in phase correlated, in the 
temporal domain, with the bolometer ones. This fact allows 
us to efficiently remove this parasitic noise via a decorrelation 
analysis in the Fourier plane using as templates the housekeep- 
ing data for the focal plane thermometers and the blind bolome- 
ter. 

Assuming a linear model we can write the Fourier trans- 
form of the parasitic signal for the bolometer / as B'""{v) = 
Fi(v) X T(v) where T represent the templates described above 
and Fi(v) is a frequency dependent correlation coefficient. As 
the parasitic signal is both well localized in time and in fre- 
quency, we perform a regression analysis in the Fourier plane 
using the windowed Fourier transform over time intervals of 
size L in order to compute the correlation coefficient F,(v). 
To maximize the efficiency of the algorithm and to limit the 
subtraction of sky signal which may be accidentally correlated 
with the templates we apply the decorrelation analysis only in 
a few predefined frequency intervals where the parasitic signal 
dominates. Further to obtain a robust estimate of the correla- 
tion coefficient we average it over time intervals so that for 
each of the frequency intervals it reads 



i:f-o'B*(v)F(y) 



Fi(y) 



(11) 



where / represents the Fourier transform of function /. 

The estimate of the high frequency parasitic noise is then 
subtracted from the data in the Fourier plane. We repeat this 
analysis for each of the available templates. The choice of the 



parameters of the method, i.e L, N, results from a trade-off : 
we want to minimize the parasitic power in the TOl spectrum, 
but we also want that the CMB signal remained unaffected by 
the process. Simulations leads to a best choice of parameters 
L - 32768 samples and = 32 to get a CMB power reduced 
by less than 1%. 

In Fig. 23 we plot the power spectra for the bolometers 
143K03, 217K06, 353K01 and 545K01 before and after sub- 
traction of the high frequency parasitic signals. We observe that 
the parasitic signal is significantly reduced in particular at the 
lowest frequencies. However for most of the bolometers around 
40 Hz the processing although efficient in removing the system- 
atic effect is not satisfactory (residuals are much larger than for 
lower frequencies) and therefore in the following the Archeops 
TOls are high pass filtered at 38 Hz before any scientific anal- 
ysis. We also observe that, in some of the bolometers very well 
localized peaks in the spectrum are not removed in our process. 
These peaks are manually characterized and the signal at their 
frequencies is removed in the Fourier domain for all Archeops 
bolometers. 



8. Data quality checks and noise properties 

As discussed in the previous section the Archeops data are af- 
fected by parasitic microphonic noise at high frequencies and 
by other unidentified systematics in the whole frequency range. 
None of them are neither stationary nor Gaussian and their 
contributions vary significantly among bolometers. In this sec- 
tion we briefly describe how we performed the selection of the 
best bolometers, in terms of level of residual systematics and 
of noise properties, which are used for the construction of the 
Archeops sky maps presented in Sect. 1 101 



8.1. Time-frequency analysis of the Archeops data 

We have performed a time-frequency analysis of the Archeops 
data using the Discrete Wavelet Packet Transform (DWPT) as 
described in |Macias-Perez & Bourrachot, (2006)| The top left 
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Fig. 24. Left column: from top to bottom, time-frequency representation of the 217T04 bolometer data and of the expected 
Galactic signal for it. Right column; same for the 217K04 bolometer. 



Fig. 25. Top: Average wavelet power spectrum in //V of the 
TOD of the 217K04 Archeops bolometer as a function of 
frequency computed from its DWPT Bottom: Time evolution 
of the power spectrum of the TOD of the 217K04 Archeops 
bolometer computed from its DWPT. 

panel of figure |23 shows the power distribution in the time- 
frequency plane for the data of the 217T04 bolometer. At low 
frequency and because of the particular scanning strategy of 
Archeops, we can observe spikes and some broader structures 
which correspond mainly to Galactic signal. This is clearly 
shown by the bottom plot where we trace the expected Galactic 
signal for this bolometer. Also at low frequency we distin- 
guish the atmospheric residuals and the l//-like noise on the 
bolometer. We can further observe from 20 to 40 Hz structures 
which vary on time being particularly strong around 21h00 UT 
time. These structures are residuals, after subtraction, of the 
high frequency systematics described before. Notice that the 
level of systematics is significantly larger than the high fre- 
quency noise making this bolometer useless for scientific pur- 
poses. 

On the top right panel of figure |23 we present the time- 
frequency analysis of the 217K04 bolometer. As above, we can 
clearly distinguish at low frequencies the Galaxy contribution 
whose template is represented in the bottom right panel. We 
observe for this bolometer no intermediate frequency struc- 
tures as for the previous bolometer. The systematic contribu- 
tion has significantly reduced by the high frequency subtrac- 
tion and the residuals are well below the high frequency noise. 
This simple but efficient analysis allows us to exclude from 
the further processing those bolometers which present either 
strong or highly time variable systematic residuals. Typically 
the Trapani bolometers presented a large amount of contam- 
ination and were excluded. We also realized that the noisier 
bolometers were in general more affected by residual system- 
atics. 

8.2. Long term non-stationarity of the noise 

On the top panel of Fig.|53] at very high frequency we can dis- 
tinguish a slow decreasing of the noise power with time which 
is the main cause of non stationarity in the Archeops data. 
This is due, as discussed in Sect. 17.11 to the sudden increase 
of temperature of the Archeops cryostat when taking off and 
to its slow cooling down to the nominal value of about 95 mK 
during the flight. To account for this non stationarity we have 
modeled the Archeops data as a time modulated-stationary 
wavelet process (Macias-P erez & Bourrac hot, (2006) I. These 
processes correspond to a continuous generalization of piece- 
wise stationary processes. They are described by two main vari- 
ables : the mean power spectrum of the data and a time vary- 



ing function, cr(t), which account for the time variations of the 
power spectrum. 

Figure |25l shows on the top and bottom panels the time av- 
eraged wavelet spectrum and the cr(f) function respectively for 
the noise of bolometer 21 7K04. On the top plot we can clearly 
distinguish the increase of noise power with frequency which 
is due to the deconvolution from bolometer time response. On 
the cr(t) function, there are two distinct regimes corresponding 
first to a fast cooling of the cryostat in the first two hours and 
then to a roughly constant temperature of the focal plane dur- 
ing the last 10 hours of flight. Notice that our estimate of cr(t) is 
noisy and therefore for further processing we generally smooth 
it. We can also observe variations on cr{t) above the noise limit 
for the second regime, that, although small with respect to those 
on the first two hours of flight, need to be taken into account in 
any further processing. 

8.3. Gaussianity of the noise 

Up to now, we have only considered the power spectrum evo- 
lution to define the level of stationarity of the data. To be 
complete in our analysis we have first to characterize the 
Gaussianity of the noise distribution function and second to 
check its evolution with time. In practice and to reduce the un- 
certainties in the analysis it is more convenient to proceed the 
other way around. 

We have implemented a Kolmogorov-Smirnov test in the 
Fourier domain to check the evolution with time of the noise 
distribution function of each of the bolometers. Notice that the 
intrinsic bolometer noise can be considered Gaussian to a very 
good approximation and therefore any changes on the distribu- 
tion function of the noise will indicate the presence of signif- 
icant residuals from systematics. We work in the Fourier do- 
main both to speed up the calculations and to isolate the noise, 
which dominates at intermediate and high frequencies, from 
other contributions like the Galactic and atmospheric signals. 
We have performed the test in consecutive time intervals of 
about 7 minutes which we compare two by two. We divided the 
frequency space in bins of 1 Hz. The test is considered to fail 
when the probability of having a greater Kolmogorov-Smirnov 
statistics under the equal distribution hypothesis is less than 
1%. 

The results of the test for the KS3 bolometers 143K03, 
217K06, 353K06 and 545K01 are shown on Fig. IH The 
white points represent failing intervals in the time-frequency 
domain. Between 16h00 and 27h00 UT, the test is successful 
except for the lowest frequency bins (f<10 Hz), in which the 
Galactic and atmospheric signals, which are neither Gaussian 
nor stationary, dominate. We notice that elsewhere the number 
of points where the test fails is no more than 1 % of the total, 
that is, exactly what is expected under the hypothesis of no 
time evolution of the distribution function. Moreover, they do 
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not exhibit any clustering. Around 28h00 UT, the sunrise on 
the gondola induces strong time evolution due to the heating 
of the 10 K stage. 



Fig. 26. From top to bottom, results of the Kolmogorov- 
Smirnov test on the bolometers 143K03, 217K06, 353K06 and 
545K01 respectively. The white polygons correspond to inter- 
vals in the time frequency plane where the test is considered to 
fail (see text for details). 



In the range from 16h00 to 27h00 UT, as the properties of 
the distribution function does not vary, we can globally study 
the Gaussianity of the noise. For this purpose we have imple- 
mented a simple test in the Fourier domain by fitting a Gaussian 
to the histogram of the coefficients of the Fourier decomposi- 
tion of the time ordered data which were binned in 1 Hz inter- 
vals as above and computing the value for the fit. We con- 
sider the test to fail if for Gaussian distributed data the probabil- 
ity of having a greater reduced that the one measured is sig- 
nificantly below 5 %. Figurel^shows the results of the test for 
the KS3 bolometers, 143K03, 217K06, 353K06 and 545K01. 
We trace the reduced as a function of the frequency bin. In 
dark and light blue we overplot the values for which the 
probability of having a larger one considering a Gaussian dis- 
tribution are 95 % and 5 % respectively. The reduced c/i/2 mea- 
sured are almost everywhere below or around the 5 % limit, ex- 
cept in the first frequency bin where the Galactic signal, highly 
non Gaussian, dominates. We can therefore consider that the 
data are in a first approximation compatible with Gaussianity. 

8.4. A/o/se power spectrum estimation and simulations 

We have performed noise simulations, which we call con- 
strained realizations of noise, to fill the gaps in the data. In 
this case we use a simple algorithm. First, we reconstruct the 
low frequency noise contribution via interpolation within the 
gap using an irregularly sampled Fourier series. Finally, we 
compute the noise power spectrum locally (in time intervals 
of about 1 hour around the gap) at high frequency and then 
we produce a random realization of this spectrum. Notice that 
we are only interested in keeping the global spectral properties 
of the data. Moreover, the gaps are in general very small in 
time compared to the piece of the data used for estimating the 
power spectrum and therefore this simple approach is accurate 
enough. 

For the first estimate of the CMB angular power spectrum 
with the Archeops data (Benoitet al. (2003a) i we needed an 
accurate estimate of the noise angular power spectrum. For this 
purpose, we have estimated the time power spectrum of the 
noise for each of the Archeops bolometers using the algorithm 
described in [Amblard & Hamilton (2 004) This algorithm re- 
lies on the iterative reconstruction of the noise by subtracting 
in the TOD an estimate of the sky signal. The latter is obtained 



from a coadded map which at each iteration is improved by tak- 
ing into account the noise contribution. From the reconstructed 
noise timeline we can then obtain for each bolometer both the 
average noise power spectrum and the cr(f) function described 
in Sect. 18.21 Since the noise in the Archeops data can be con- 
sidered as Gaussian distributed (Sect. l83t these two quantities 
are enough to simulate noise timelines using the algorithm pre- 
sented in Macias-Perez & B ourrachot, (2006)| The fake time- 
lines can then be projected into maps for the estimate of the 
noise angular power spectrum. In general a few hundreds sim- 
ulations are needed to obtain a reasonable estimate. 

8.5. The best bolometers 

From the above analysis of the Archeops noise we have chosen 
for each frequency band those bolometers for which the resid- 
ual systematics are well below the noise level. For the CMB 
channels at 143 GHz and 217 GHz the selected bolometers 
(143K03, 143K04, 143K05, 143K07, 217K04 and 217K06) 
correspond to the more sensitive ones going from 94 to 200 
HK .s'^^. At 353 GHz the bolometers were all of similar quality 
and we keep all of them. At 545 GHz only one bolometer was 
available, 545K01. 

9. Calibration 

We describe in this section the global absolute calibration and 
intercalibration of the Archeops data. The former is performed 
using three different types of calibrator : the CMB dipole, the 
Galaxy and the planets Jupiter and Saturn. At low frequencies 
(143 and 217 GHz) the dipole is the best absolute calibrator. 
At higher frequencies we need to use the Galaxy because the 
dipole signal is too faint with respect to the noise and system- 
atics. 

9.1. CMB dipole 

At low frequencies, the CMB dipole is a very good absolute 
calibrator jPiat et al. (2003)1 [Cappellini et al. (2003)t and 
therefore it constitutes the primary absolute calibration of the 
Archeops data. Here we use the total dipole which is the sum 
of the solar dipole (constant in time) and the Earth induced 
dipole (with annual variations due to the Earth change of 
velocity) computed at the time of flight. Indeed, the dipole 
calibration has the following advantages: 1) the dipole is 
spread over the all sky and thus it is always present whatever 
the pointing, 2) it is much brighter (typically a factor of 100) 
than the CMB anisotropics, but still faint enough so that 
non linearity corrections are usually not needed, and finally 
3) it has the same electromagnetic spectrum as the CMB 
anisotropics so that no color corrections need to be applied. 
The only drawback is that we must assume a constant response 
for the instrument throughout a wide range of angular scales, 
i.e. an extrapolation from ^ = 1 to 1000. The dipole being 
an extended source we need to account for the beam and the 
spectral transmission of each of the detector to generate the 
point source calibration. 
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Fig. 27. From left to right and up to bottom, maximum reduced (97 d.o.f) of the Gaussian fit of the histogram of the Fourier 
coefficients of the Archeops data for the 143K03, 217K06, 353K06 and 545K01 bolometers respectively. In dark and light blue 
we overplot the values for which the probability of having a larger one considering a Gaussian distribution are 95 % and 5 % 
respectively, (see text for details on the analysis) 



In the Archeops case, the dipole signal is expected to 
contribute only to the fundamental rotation frequency /spin. 
However, the drift of the rotation axis during the flight leads to 
a broadening of the /sp„, dipole contribution which needs to be 
taken into account. Furthermore, as discussed in Sect. l7.2l other 
signals show up at the spin frequency including both Galactic 
and atmospheric signals. To account for these, we compute the 
total dipole calibration coeflicients from the correlation analy- 
sis described in Sect. 17.21 We actually produced templates of 
the Galactic emission using extrapolation of the dust emission 
to the Archeops frequency ( Finkbeiner et al (1999)1 and tem- 
plates for the atmospheric emission using the housekeeping 
data from the altitude and elevation of the balloon. The tem- 
plate for the total dipole in time domain was obtained by depro- 
jecting, using the Archeops pointing, a simulated CMB dipole 
map for which the dipole amplitude and direction were taken 
from the WMAP results (Bennett et al. (2003)1. The left panel 
of Fig. shows the simulated total dipole map in Galactic 
coordinates for the Archeops coverage and centered in the 
Galactic anticenter On the right panel we plot the Archeops 
reconstructed CMB dipole for the 143K03 bolometer. We ob- 
serve a very good agreement between the two maps but for 
some stripes when crossing the Galactic plane and on the North 
at high Galactic latitude. 

Indeed, the CMB dipole is detected with a signal to noise 
better than 500 for 12 hours of data on the 143 and 217 GHz 
channels. The errors on the dipole calibration coefficients come 
mostly from systematic effects. We produced several versions 
of the calibration coefficients by changing the templates used 
in the fit. We noticed that the result is very stable with respect 
to even the most extreme cases. From those tests, we can de- 
duce that the overall uncertainty (1 cr) is 4 % (resp. 8 %) for 
the 143 (resp. 217) GHz photometric pixels. The larger uncer- 
tainty for the 217 GHz channel reflects the spectrum of the 
atmospheric contamination. Notice that the dipole calibration 
is performed at the early stages of the analysis before remov- 
ing the spin frequency systematics (see section 0. When ap- 
plied to the 353 GHz channel, there is a residual contaminant 
at the same level as the dipole. Considering the 545 GHz chan- 
nel in the analysis can help in this case and allows finding the 
dipole calibration coefficient at 353 GHz in agreement with the 
Galactic calibration coefficients within 20 % as described be- 
low. 

9.2. The Galaxy 

At high frequencies 353 and 545 GHz, the Galactic emission 
has to be used to calibrate the Archeops data. The best data, 
in term of spectral coverage and absolute calibration accuracy, 
are the FIRAS spectra (Mather et al. 1990) . FIRAS, on board 



the COBE satellite, is a scanning Michelson interferometer 
that has provided low-resolution spectra in a low (2 - 20 
cm"') and high (20 - 100 cm"') frequency band, for 98.7% 
of the sky. For individual pixels, the signal-to-noise ratio of 
the FIRAS spectra is about 1 at high Galactic latitude and 
~50 in the Galactic plane. The FIRAS maps used here were 
obtained by fitting each FIRAS spectrum using a modified 
Black Body with a \^ emissivity law and extrapolating the fit 
at Archeops frequencies. Since we are searching for the best 
representation of the data and not for physical dust parameters, 
we restricted the fit to the frequency range of interest (this 
avoids the need for a second dust component, of the type 
proposed by Finkbeiner e t al (1999) I. The FIRAS brightness 
maps are then converted to photometric maps with the flux 
convention of constant v/y. To be compared to the FIRAS maps 
and to maintain the best possible photometric integrity, the 
Archeops data are convolved by the FIRAS beam and put into 
the FIRAS Quadrilateralized Spherical Cube (QSC). At this 
stage, we have Archeops maps at the FIRAS resolution that 
can be directly compared to the FIRAS maps at the Archeops 
frequencies to derive the Galactic calibration factors. 

Due to its high signal-to-noise ratio and its extension, the 
Galactic plane is the best place to derive the calibration factors. 
We compute Galactic latitude profiles of both maps (\b\ < 30°) 
at selected longitudes and perform a best straight-line fit to 
the Archeops-FIRAS profile correlation from which we derive 
the calibration factor and its error bar An example of this fit 
is shown on Fig.l29lwhere we observe a linear correlation be- 
tween the Archeops and FIRAS profile intensities. The whole 
calibration process gives calibration factors with statistical 
errors of about 6 %. Details on the whole calibration process 
are given can be found in |Lagache & Douspis (2006)| 



Fig. 29. FIRAS/Archeops Galactic profiles correlation on the 
Galactic plane (bolometer 353K01 at 353 GHz). Fitting a 
straight-line gives the calibration factor 



Our procedure has been extensively tested. We have 
applied the calibration scheme to the comparison between the 
FIRAS and DIRBE data at 140 and 240 yum and found results 
in very good agreement with those of Fixse n et al (1997)| 
We also tested our procedure on the Finkbein er et al (1999)| 
maps, although these maps exclude the Galactic plane below 
\b\-l°. We obtain gain differences which are less than 6% 
across the Archeops frequencies. We also used Archeops maps 
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Fig. 28. From left to right, simulated and reconstructed CMB dipole maps for the Archeops 143K03 bolometer centered at the 
Galactic anticenter. To reconstruct both dipole maps the timelines have been band pass filtered. This introduces discontinuities 
on the maps due to the scanning strategy. 



obtained using different low-frequency filters to test the effect 
of the filtering on the calibration. We found that it modifies the 
calibration factors by only ~3%. 

Although the calibration on the Galaxy is less accurate 
for the low-frequency than for the high-frequency channels of 
Archeops, we compare the CMB dipole and Galaxy calibration 
factors on Fig.|30| We observe a good agreement between both 
calibrations, with a mean ratio of the calibration factors lower 
than 1.05 and all of them compatible with 1 within the Icr error 
bars. This clearly demonstrates the robustness of the methods 
used for the extended emission calibration and the consistency 
of the data reduction from TOls (dipole calibration) to maps 
(Galactic calibration). 



Fig. 30. Comparison of the Galactic and Dipole calibration 
factors (in mK///V). The error bars (both statistical and sys- 
tematic) are of about 4% and 8% for calibration on the dipole 
at 143 and 217 GHz respectively and from 6% to 15% for the 
calibration on the Galaxy. 



9.3. Jupiter calibration 

We have also performed a point source calibration using the 
measurements at the two independent crossings of the plan- 
ets Jupiter and Saturn (5 times fainter than Jupiter). For this 
purpose we proceed as follows: 1) We obtain the beam pat- 
tern shape from the Jupiter measurements after deconvolution 
from the bolometer time response 2) We compute the flux of the 
sources (in fiW) using pixel photometry up to a radius of 40 ar- 
cmin starting from the center of the beam pattern. 3) Finally, 
we compare the fluxes to a model of the temperature emission 
of the sources ( [Moreno (199 8) l that reproduces radio observa- 
tions of Jupiter and Saturn from 20 to 900 GHz. Taking into 
account uncertainties on the model, absolute errors on point 
source calibration with Jupiter are estimated to be ~12%. 

The ratio between Jupiter and Saturn flux measurements at 
a given frequency does not depend on the absolute instrument 
calibration. We find ratios of 0.97 ± 0.016 and 1.02 + 0.020 at 
143 and 217 GHz respectively which are not compatible with 
flie one measured by [Goldin et al. (1997)| (0.833 + 0.012) at 
similar frequency (171 GHz). Due to the large brightness of 
Jupiter (about 1 K^y equivalent brightness) we expect some 
non-linearity on the bolometer response which could be the 
cause of this difference. This problem in addition to the uncer- 
tainties on the knowledge of the beam pattern and in particular 
of the far side lobes could also explain the small discrepan- 



cies between the Jupiter calibration and the dipole calibration 
as shown on Fig. 13 II 

Fig. 31. Comparison of the point source (Jupiter and Saturn) 
and dipole calibration factors. The error bars are about 4% and 
8% for calibration on the dipole at 143 and 217 GHz respec- 
tively and about 12% for the calibration on the sources (essen- 
tially due to the uncertainty of the thermal emission model). 



9.4. Intercalibration 



Bolometer 


relative calibration (fiV/fiW) 


error (stat.) 


143K01 


1.16 


2.1% 


143K03 


1 




143K04 


1.63 


2.3% 


143K05 


1.15 


1.8% 


143K07 


1.39 


1.9% 


143T01 


2.06 


2.3% 


217K01 


1.80 


1.1% 


217K02 


1.72 


0.8% 


217K03 


9.09 


2.5% 


217K04 


0.903 


0.7% 


217K05 


3.21 


1.2% 


217K06 


1 




217T04 


1.50 


1.1% 


217T06 


1.02 


1.02% 


353K01 


1 




353K02 


1.13 


0.68% 


353K03 


1.19 


0.70% 


353K04 


1.02 


0.70% 


353K05 


1.20 


0.79% 


353K06 


1.05 


0.78% 



Table 6. Relative calibration coefficients and their relative sta- 
tistical errors from Galactic profiles at constant Galactic longi- 
tude. 



The polarization signal for experiments like Archeops is re- 
constructed, in a first approximation, from the differences be- 
tween pairs of bolometers (see Sect. 110. 3> . Therefore the ac- 
curacy on this reconstruction is very sensitive to the relative 
calibration between bolometers. In the case of the Archeops 
353 GHz polarized channels, the absolute calibration on the 
Galaxy presented above is not accurate enough for the di- 
rect reconstruction of the polarized maps. Thus, we have im- 
plemented a relative caUbration algorithm based on the inter- 
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comparison, for all bolometers in the same channel, of Galactic 
profiles at constant Galactic longitude. This algorithm has also 
been applied to the unpolarized channels at 143, 217 GHz as a 
cross check of the absolute calibration analysis. 

Assuming bolometers per channel we can measure, N 
Galactic profiles as a function of the Galactic latitude, b, 

Sj(b) - ajs(b) + nj(b) 

where aj and nj(b) are the intercalibration coefficient and 
noise contribution for the bolometer j respectively, and 
is the true Galactic profile at the frequency of interest. Our 
algorithm, which is based on a minimization constrained 
via Lagrange's multipliers, estimate simultaneously the 
intercalibration coefficients and the mean Galactic profile. 
This algorithm is described in details on appendix |X| We 
have performed robustness test to validate the algorithm. In 
particular we have checked that the results are not affected 
by the choice of constraint and that the analytical error bars 
obtained for the intercalibration coefficients are reliable. 



Fig. 32. From left to right, Archeops Galactic profiles in fiV 
at constant Galactic longitude respectively before and after the 
intercalibration for the 353 GHz bolometers. 



The algorithm was applied to the Archeops data which 
were previously processed as described in the previous sec- 
tions. The profiles were obtained by averaging the signal sam- 
ples within latitude bins. The errors are computed assuming 
that the noise is white, uncorrelated between bolometers and 
stationary. The noise of each sample is estimated from either 
data outside the Galactic plane (\b\ > 25°) or high-pass filtered 
data above 10 Hz (where Galactic signal is negligible). Results 
are very similar (differences between the two methods are less 
than 3%, inducing negligible difference in intercalibration co- 
efficients) in the two cases. Finally, care has been taken on 
keeping the same sky coverage for all bolometers when com- 
puting the Galactic profiles. 

For the polarized channel the presence of strongly polarized 
regions on the sky may affect the computation of the intercal- 
ibration coefficients. To avoid this, we proceeded in two steps. 
First, we compute the relative calibration coefficients and mean 
Galactic profile on the full common sky area for the six polar- 
ized bolometers. Using these coefficients we build polarization 
maps and label strong polarized areas which are then excluded 
from the analysis. We then build again the Galactic profiles and 
recompute the relative calibration coefficients. After two iter- 
ations, we observed that the estimates of the intercalibration 
coefficients are stable. We show on Fig. |^ the results of this 
analysis for the 353 GHz channel. From left to right we plot 
the Galactic profiles before and after intercalibration respec- 
tively. We observe that the intercalibration is achieved with a 
high degree of precision. 

Final results on the intercalibration coefficients for the 
Archeops data are given in Tab. |6| for all the bolometers at 



143 GHz, 217 GHz (unpolarized) and 353 GHz (polarized). 
The Icr statistical errors on the relative calibration coefficients 
are at most 2.5% and always below 1% for the polarized 
channel. Error bars are smaller at high frequency because the 
Galaxy signal is stronger. 

9.5. Overall sensitivity 

Using the observed response and white noise level of the 
bolometers, we can now estimate the global efficiency of the 
instrument during the KS3 flight. Table^lgives the instrument 
sensitivities in various units at the timeline level. This places 
Archeops within the best range of instrumental instantaneous 
sensitivities for millimetre continuum measurements. We also 
compute the averaged sensitivity, quoted in Tab. |8j by using a 
total integration time of 12 hours and a total sky coverage of 
30 %. In this case Archeops is not within the best experiments 
because the sensitivity per pixel is diluted by the large covering 
area needed for reconstructing the large angular scales on the 
sky. 

10. The Archeops sky maps 

In this section we finally describe how we obtain submillimetre 
sky maps from the Archeops timelines and pointing informa- 
tion. First and prior to projection we remove low frequency 
drifts on the Archeops timelines via a destriping algorithm. 
Then, these timelines are processed in three different ways to 
produce after projection: CMB, Galactic intensity and polar- 
ization maps. 

10.1. Destriping 

Even after the subtraction of the very high and low frequency 
identified systematic effects from the data, we can observe 
residual stripes on the Archeops simple coadded maps. These 
are mainly due to low frequency drifts in the data coming 
from atmospheric residuals and other thermal backgrounds. 
Notice that there is no electronic 1 // component because of 
the AC bolometer modulation scheme (see Sect, \4-.2i . To re- 
move those drifts we have implemented a destriping algorithm 
l |Bourrachot (20 04)^ making the assumption that the scanning 
direction is generally not related to the orientation of the struc- 
tures on the sky. To destripe we compute a low frequency base- 
line in the timelines by minimizing the ratio between the rms 
in the cross-scan and in the in-scan directions directly from 
the time ordered data (no reconstruction of maps is needed). 
To represent the baseline we use a basis of localized functions 
Uk{t) where 

U{t) = ./„c(f)exp(-5^) 

Udt) = Uit-t,) (12) 
til - kA. 

These functions are regularly sampled and contain only fre- 
quencies lower than 1/(2A). 

The minimization is performed outside the Galactic plane 
over boxes with sizes related to A (typically a few tens of square 
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Freq.(GHz) 


Nboi 


A/Krj/Hz"- 


MJy/sr/Hz'/^ 


A/Kcmb/Hz"^ 


IOVHz"^ 


/jKcmb.s"- 


143 


6 


50 


0.031 


87 


23 


61 


217 


7 


39 


0.056 


127 


286 


90 


353 


6 


82 


0.315 


1156 


132 


817 


545 


1 


77 


0.702 


9028 


495 


6384 



Table 7. The Archeops KS3 in-flight timeline sensitivities per channel. The best bolometers on each frequency channel are op- 
timally combined {Nhoi) to obtain a sensitivity at the instrument level. From left to right the units correspond to a Rayleigh-Jeans 
spectrum, then a constant v/,, spectrum, and a CMB spectrum. Sensitivity to the SZ effect is measured with the dimensionless y 
Compton parameter To convert to 1 cr and to one second integration, we can simply divide by V2: see an example in the last 
column. 



Freq.(GHz) 


Ntoi 




MJy/sr/Hz 




lO^v 


Jy 


143 


6 


57 


0.036 


98 


27 


1.2 


217 


7 


44 


0.064 


144 


325 


2.2 


353 


6 


94 


0.358 


1312 


150 


12.1 


545 


1 


87 


0.797 


10251 


562 


27.0 



Table 8. The Archeops KS3 in-flight map sensitivities per channel. The best bolometers on each frequency channel are optimally 
combined {Nhoi) to obtain a sensitivity at the map level. A square pixel of 20 arcminutes is taken to compute the average 1 cr 
noise. KS3 flight roughly covered 30 % of the total sky, which represent 1 10,000 pixels. A bolometer has observed a pixel on the 
map during an average time of 0.4 seconds. 



Fig. 33. Left: Power spectrum of the time ordered data of the bolometer 545K01 before (black) and after (red) destriping. Right: 
Zoom-in of the left plot at first multiples of the spinning frequency. 



degrees). The cut on the Galactic plane is obtained from a 
Galactic mask derived from Galactic template maps computed 
at the Archeops frequencies (using Finkbeiner et al (1999) ). 
The algorithm is applied in steps of decreasing A (4000, 1000, 
500 and 300) to focus on different frequency ranges. In any 
case, data at frequencies above 1 Hz are not affected by this 
method. Prior to this process we generally apply a classical de- 
striping algorithm based on the minimization of the variance 
per pixel in the maps to produce a first approximation of the 
baseline and in particular of the lower frequency components 
(below 0.7 Hz). 

We have applied the destriping algorithm to simulated 
Archeops data at 217 and 545 GHz including correlated noise 
at low frequencies and Galactic and CMB emissions. No bias 
has been observed in the estimation of the Galactic signal. For 
the CMB the power is reduced by at most 5 %. The full de- 
striping transfer function in the multipole space is presented in 
[Tristram et al. (2005b)| 

An example of the application of the destriping procedure 
to the Archeops data is shown on Fig. |33] On the left panel, 
we represent the power spectrum of the 545K01 bolometer 
data before (black) and after (red) application of the destriping 
method. We observe that the noise is reduced significantly at 
low frequencies and in particular the spectrum flattens. Further 
we can observe from the right plot that power spectrum signal 
at the few first multiples of the spinning frequency are much 
broader before destriping. We can conclude, by comparing the 
353 and 545 GHz data, that this extra structures come mainly 
from the atmospheric emission. 



10.2. Specific processing for Galactic maps 

10.2.1. Atmospheric contamination 

For producing Galactic maps from the Archeops data we need 
first to remove the residual parasitic atmospheric noise. The 
destriping algorithm described above, although very efficient 
at frequencies lower than 1 Hz, can not fully eliminate it. From 
Fig.|^we observe that the latter shows up on the power spec- 
trum of the time ordered data as residuals at the spin frequency 
multiples. This also produces two well defined structures, be- 
tween 0.9 and 1.6 Hz, in the power spectrum. This is shown 
on Fig. where we plot the power spectrum of the TOI for 
the 143K03, 217K04, 353K01 and 545K01 bolometers. Notice 
that the parasitic noise shows a common spectrum shape for all 
the Archeops bolometers and that its total intensity increases 
with the frequency of observation. 

To estimate this parasitic atmospheric noise and re- 
move it from the data we have used a modified ver- 
sion of the MCMD-SMICA component separation algorithm 
( Delabrouille et al. (2003)) which can work directly on time or- 
dered data. We have assumed a very simple linear model for 
the Archeops timelines with three main components: Galactic 



Fig. 35. Power spectrum in arbitrary units of the parasitic-like 
component for the MDMC analysis of the bolometer 353K01 
for different time intervals. 
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Fig. 34. From left to right and from top to bottom MDMC decomposition of the Archeops data at intermediate frequencies for 
the 143K03, 217K04, 353K01 and 545K01 bolometers. The black, blue and red line correspond to power spectrum in arbitrary 
units of the raw data, the parasitic-Uke and the Galactic-like contributions respectively. 



emission, atmospheric emission and Gaussian instrumental 
noise. In the time ordered data the Galactic emission is weak 
relative to the noise. Therefore, to improve the convergence of 
the algorithm we have used as inputs, apart from the Archeops 
data, fake timelines of Galactic emission extrapolated to the 
Archeops frequencies from the IRAS maps using model 8 in 
[Finkbeiner et al (1999)| Finally, to reduce the noise contribu- 
tion we have restricted our analysis to the frequency range 
from 0.03 to 2.5 Hz. The main results of this analysis for the 
143K03, 217K04, 353K01 and 545K01 bolometers are pre- 
sented in Fig. The blue and red curve are the reconstructed 
atmospheric and Galactic emissions. We observe that the re- 
constructed emission reproduces very well the two structures 
on the power spectrum and also contributes to the multiples 
of the spin frequency. As we expect the atmospheric emission 
to vary with time we have performed this analysis for differ- 
ent intervals. As shown in Fig.|35]the power spectrum of the 
atmospheric parasitic noise does not change significantly nei- 
ther in shape nor in intensity. From these results we have con- 
structed a template of the atmospheric emission which is sub- 
tracted from the data via a simple decorrelation analysis as de- 
scribed in Sect.0 

10.2.2. Galactic maps 

The final Archeops Galactic maps, presented in Fig. |36l are 
produced in the Healpix pixelization scheme by simple coaddi- 
tion of the previous processed timelines which are previously 
band-pass filtered. The low-pass filtering allows us to both re- 
move spurious high frequency noise in the data (see section^ 
and avoid aliasing on the final maps. The high-pass filtering 
keeps frequencies above 0.03 Hz and to reduce ringing we first 
mask the brightest Galactic regions and fit an irregularly sam- 
pled Fourier series truncated to the frequency of interest. The 
latter is then fully sampled and subtracted from the data. We 
produce individual maps for each of the detectors as well as 
combined maps per channel using the best available bolome- 
ters. 

From top to bottom. Fig. |36| shows the combined Galactic 
maps for the Archeops 143, 217, 353 and 545 GHz channels 
respectively. These are the first available large angular scales 
maps of the sky in this frequency range. The maps are displayed 
in antenna temperature units and in Galactic coordinates with 
the Galactic anticenter at the center of the map. The Galactic 
plane structure, including for example the Cygnus region on the 
right of the map and the Taurus region on the left, are clearly 
visible. Their intensity increases globally with frequency as ex- 
pected for Galactic dust emission. At high Galactic latitudes 
the maps at low frequencies show no contamination from the 
atmospheric emission. At 545 GHz we can observe some at- 
mospheric contamination. This was expected since the atmo- 



spheric signal is stronger at high frequencies and because we 
have only a single bolometer available. More detailed descrip- 
tion and scientific analysis of these maps will be presented in a 
forthcoming paper In particular almost all (about 100) identi- 
fied point sources are Galactic and will be presented elsewhere. 

10.3. Specific processing for polarization maps 

Whereas most of the preprocessing and noise subtraction 
is common to all channels, the 353 GHz polarized channel 
requires additional specific treatments. The direction of 
polarization of the bolometers oriented mechanically in the 
focal plane has to be checked, as well as their polarization 
efficiencies. Moreover, the map making algorithm differs from 
that of temperature maps since it has to deal with non scalar 
quantities. Details regarding the map making algorithm as well 
as the final polarized Archeops maps at 353 GHz are given 
in [Benoit et al. (2004)| Previous to map making the TOIs are 
processed as above using the MCMD-SMICA algorithm to 
remove the contamination from atmospheric emission. 

The reconstruction of the polarization directions of the 
bolometers in the focal plane was performed during ground cal- 
ibration. For this purpose, we built two wire grid polarizers of 
10 cm diameters with the same technology as the one used for 
wire chambers in high energy physics experiments. We used 
Cu/Be wires of 50 microns, spaced by 100 microns on circular 
steel frames which was expected to produce an incoming ra- 
diation at more than 98% at 850 microns (353 GHz). We built 
an alignment mechanism that could hold one or both polarizers 
on top of the entrance window of the cryostat, facing directly 
the focal plane. One of the polarizers could rotate at 1.5 rpm. 
A 13.4 Hz chopper modulated the incoming radiation of a liq- 
uid Nitrogen polystyrene box used as a 77 K black body to 
enable a lock-in detection. Once the lock-in and standard noise 
subtraction were performed, the rotating polarizer induced a 
1.5 rpm period sinusoidal response for the polarized channels, 
from which the phase provided the direction in the focal plane, 
and the off'set and amplitude the cross-polarization level. The 
positions were confirmed to be nominal (30, 120, 150, 240, 90, 
0) deg w.r.t. the scan axis up to the precision of the method 
which was estimated to be 3 deg. The cross polarization, de- 
fined as the ratio cross-Intensity/co-Intensity was found to be 
approximately 2%. 

The absolute calibration of the polarized photometers was 
performed in the same way as for the other channels. The rela- 
tive calibration was performed using the algorithm described in 
Sect. 19.41 Actually, it was originally designed for the polarized 
channels. The only caveat was the possible systematic effect 
induced by a polarized Galactic component. A two step itera- 
tion process was designed, the first of which consisted in the re- 
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Fig. 36. From top to bottom: Galactic maps in antenna temperature for the 143, 217, 353 and 545 GHz Archeops channels. They 
are displayed in Galactic coordinates with the Galactic anticenter at the center of the map. 



Fig. 37. From top to bottom, power spectrum of the Archeops 
time ordered data before (black curve) and after (red curve) 
foreground removal for the 143K03 and 217K04 bolometers 
respectively. For comparison, the bottom plot shows the power 
spectrum of the 545K01 bolometer. 



Fig. 38. Top: Angular power spectrum of the simulated CMB 
signal before (black line) and after data processing (blue line) 
for the 143 GHz. Bottom: Transfer function of the data process- 
ing pipeline for 143 GHz data. 



moval of the strong polarized regions before a final relative cal- 
ibration on the non polarized parts of the Galaxy. This process 
ensured a relative calibration of the polarized detectors better 
than 2%. More details can be found in |Benoit et al. (2004)| 

10.4. Specific processing for CMB maps. 
10.4.1. Foreground removal 

For reconstructing the CMB signal on the sky we need to re- 
move from the Archeops data the foreground contribution cor- 
responding to the Galactic dust emission and the atmospheric 
parasitic noise. Both of them have a rising electromagnetic 
spectrum with increasing frequency. Therefore they are sig- 
nificantly brighter in the Archeops high frequency channels. 
The Galactic dust emission in the millimetre and submillime- 
tre range has a grey body spectrum with an emissivity of 
the order of between 1.7 and 2 (see Fin kbeiner et al (1999)| 
[Bennett et al. (2003)| [Lagache (2003) for more details). A 
slightly steeper index increasing with frequency is observed for 
the atmospheric parasitic noise. 

We have developed a decorrelation algorithm to remove 
the foregrounds from the Archeops data. As templates of the 
foreground emission we have used the Archeops high fre- 
quency channel bolometers but also fake timelines of the ex- 
pected Galactic emission contribution to the Archeops data 
to improve the efficiency of the algorithm. These fake time- 
lines were produced in two steps. First we extrapolated the 
IRAS satellite data to the Archeops frequencies using the 
model 8 of Fink beiner et al (1999)| Second the extrapolated 
maps were deprojected into time ordered data following the 
Archeops scanning strategy. For the decorrelation analysis of 
each of the low frequency (at 143 or 217 GHz) bolometers 
we used as templates the 353 and 545 GHz time ordered data, 
fake Galactic timelines corresponding to those data and an ex- 
tra fake Galactic timeline corresponding to the decorrelated 
bolometer. 

To improve the efficiency of the decorrelation method we 
bandpass filter both the Archeops data and the fake Galactic 
timelines in the range 0.1 to 2 Hz where the atmospheric and 
Galactic emission dominate. This can be clearly seen on the 
bottom plot of Fig.|^where we represent the power spectrum 
of time ordered data for the bolometer 545K01. At low fre- 
quency we observe the Galactic and atmospheric emissions in 



the form of peaks at frequencies which are multiples of the 
spinning frequency. Between 1 and 1 .5 Hz we can distinguish 
the atmospheric parasitic structure discussed in Sect. 110.21 
Above 1.6 Hz the instrumental noise dominates. The correla- 
tion coefficients are computed via a simple regression analysis 
using the bandpass filtered data. A linear combination of the 
templates previously smoothed and multiplied by the correla- 
tion coefficients is removed from the data. 

The first two upper plots of Fig. |^ show the results of 
the decorrelation analysis for the Archeops bolometers 143K03 
and 217K04 respectively. We plot the power spectrum of time 
ordered data before (black curve) and after (red curve) decor- 
relation. For the 143K03 bolometer we can observe that the 
peaks in the spectrum are completely removed by the decor- 
relation analysis but at frequencies lower than 0.2 Hz. The 
same is found for the 217K04 bolometer Further, we see that 
the atmospheric structures between 1 and 1.5 Hz are also re- 
moved. The residual Galactic emission at frequencies lower 
than 2 Hz increases dramatically with decreasing frequency. 
This seriously limits the size of the largest angular scale for 
which the CMB angular power spectrum can be reconstructed 
using the Archeops data. Although the algorithm is very effi- 
cient residual atmospheric and Galactic emission are expected 
in the final Archeops CMB maps. A more detailed discus- 
sion of these two issues is given in [Tristram et al. (2005b)| and 
IPatanchon et al. (2005)| 

10.4.2. The pipeline transfer function 

For an accurate determination of the CMB power spectrum 
with the Archeops data we have to correct from the bias in- 
troduced in the signal by the data processing. For this we have 
estimated the Archeops pipeline transfer function in multipole 
space. Here we consider the full data processing but the de- 
striping for which the transfer function was discussed above. 

For each of the Archeops bolometer at 143 and 217 GHz 
we have obtained fake Archeops CMB timelines. These were 
produced from the deprojection, using the Archeops pointing 
solution, of the same simulated CMB map. These CMB time- 
lines have been converted into voltage units using the standard 
calibration coefficients for each bolometer and then filtered out 
with the low pass prefilter (see Sect. 14. U . Further, we have 
added, to each original bolometer time ordered data, the cor- 
responding fake CMB timeline and then saved them into files 
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the same way the true ones are. Finally, we have passed each 
combined timeline through the full Archeops data pipeline but 
through the destriping. The effect of the data processing on the 
CMB simulated signal can be easily obtained. First we sub- 
tract from the combined processed timeline the corresponding 
equally processed Archeops data. Finally, we reproject the dif- 
ference into a map. Comparing the CMB angular power spec- 
trum of the simulated Archeops CMB signal before and after 
data processing we obtain the pipeline transfer function. 

The top panel of Fig. |38l shows the angular power spec- 
trum for the simulated CMB data before (black line) and after 
(blue line) data processing for one of the 143 GHz Archeops 
bolometers. Dividing the one by the other we can estimate the 
Archeops pipeline transfer function which is shown in the bot- 
tom panel plot. We observe that the changes induced in the 
CMB signal by the data processing pipeline are smaller than 
1%. Similar results are obtained for the other bolometers at 
143 and 217 GHz. Therefore, there is no need to account for 
the pipeline transfer function when estimating the CMB angu- 
lar power spectrum with Archeops (see |Benoit et al. (2003a)1 
[Tristram et al. (2005b)> . 

10.4.3. Archeops CMB maps 

The Archeops CMB maps were obtained by projection and 
bandpass filtering of the foreground cleaned timeUnes using 
the Mirage optimal map making code ( |Yvon & Mayet (2005 )^ . 
The data were low pass filtered at 30 Hz to remove spurious 
high frequency noise and high pass filtered at 0. 1 Hz to reduce 
the contribution from residual atmospheric and Galactic emis- 
sions. We have produced both individual maps for each of the 
selected best bolometers at 143 and 217 GHz, and combined 
naturally weighted maps for each of the low frequency chan- 
nels. 

Figure |39l shows from top to bottom, the Archeops com- 
bined CMB maps at 143 and 217 GHz. These maps are in longi- 
tude rotated Galactic coordinates with the anticenter at the cen- 
ter of the maps. We represent them in CMB temperature units. 
Notice that close to the Galactic center and in particular near by 
the Cygnus region (right of the map) we can observe residuals 
from the Galactic emission. However at high Galactic latitude 
neither atmospheric nor Galactic residuals are observed. For 
CMB analysis ( |Benoit et al . (2003a) Tristram et al. (2005b) i 
we use a Galactic mask to exclude the observed contami- 
nated regions. This mask is overplotted in blue on the fig- 
ure. A detailed analysis of the properties of these maps and 
how they compare to those from the WMAP satellite is pre- 
sented in a forthcoming paper ( Patanchon ef a/. (2005)| l. Other 
than CMB studies, these maps were used in combination 
with the WMAP data (Bennett et al . (2003)t to study statis- 
tically the Sunyaev-Zeldovich effect in clusters of galaxies 
( Hernandez -Monteagudo et al. (2006)^ . 

11. Conclusions 

We have presented in this paper the full processing of the 
Archeops data, from the raw telemetry to the final sky maps. 
Despite intense preparatory work, most methods and proce- 



dures discussed here were developed and implemented after 
the acquisition of the real flight data. This was mainly due to 
the requirement of a large sky coverage in a short total integra- 
tion time (about 24 hours) which imposed a large-circles like 
scanning strategy with small redundancy and therefore making 
systematic eff'ects diflicult to handle. Typically the data were 
contaminated by the large scale fluctuations of the atmospheric 
emission and by the Galactic foreground emission. Because of 
these difficulties we were forced to apply different processing 
techniques to the data for each of the main scientific goals, 
1) estimation of the CMB temperature anisotropics power 
spectrum, 2) study of the Galactic diffuse emission and 3) esti- 
mation of the polarized submillimetre emission of the Galaxy. 
In particular, the destriping and filtering techniques previous to 
projection on the sky were different for each pipeline leading 
to different output maps. A common general destriping, based 
on the assumption that the structures do not have preferred 
directions on the sky, was applied to data in all pipelines. For 
CMB maps, the low frequency channels were further decor- 
related from a mixture of the high frequency data, which are 
dominated by atmospheric and Galactic signal. For Galactic 
maps, the atmospheric component was subtracted using a com- 
ponent separation method on the timelines. For the polarization 
pipeline simultaneous time and frequency filtering was applied. 

The processing and the instrumental setup were improved 
between successive flights going from the Trapani test flight to 
the latest Kiruna one. For example, after analysis of the data 
of the first two flights, we were able to reduce significantly the 
high frequency noise excess in the data by moving the spinning 
pivot motor higher up in the flight chain to dampen mechanical 
vibrations. A thermal dependency of the signal with the 10 K 
thermal stage was completely removed for the last flight after 
a complete clean-up of the corrugated back-to-back horns in 
each of the Archeops photometric pixels. By contrast, because 
of the wearing off of the instrument from flight to flight and the 
lack of time for a complete ground calibration, the instrument 
was not launched in a fully optimal configuration in the latest 
Kiruna flight. In particular, we are aware that as a consequence 
of the accidental landing during the penultimate campaign, 
the telescope was slightly out of focus producing optical 
beams larger and more elongated than expected. To correct 
for this asymmetry new specific processing techniques were 
developed. 

The processing of the Archeops data could not be per- 
formed in a single linear pipeline. We needed extra pipelines 
to reconstruct the pointing information, compute basic in- 
strumental and observational parameters. For example with 
respect to calibration, we had to design dedicated pipelines 
for each of the method used: dipole reconstruction. Galaxy 
dust emission, intercalibration, planet calibration. Equally, 
we have worked in independent pipelines to characterize the 
instrumental response but also the data processing in order to 
correct the bias introduced by those effects in the computation 
of the CMB temperature angular power spectrum and of 
the polarized dust power spectra. For this purpose we have 
computed the transfer functions in multipole space for the 
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Fig. 39. From top to bottom, combined Archeops CMB maps for the 143 and 217 GHz channels. In the Galactic plane region the 
residual galactic emission is still visible but clearly disappears at high galactic latitudes where the CMB studies are performed. 



beam smoothing, the pipeline processing and the destriping. 
Further, we have produced very specific pipelines to remove 
foreground contributions on the data for each of the scientific 
goals. Many checks were performed on the subpipeline levels 
using simulated data. Complete end-to-end tests were difficult 
to achieve because of the complexity of the problem and the 
fact that the pipelines could not be gathered in a single one. 

Archeops provides the first submillimetre maps of the sky 
with large sky coverage, of the order of 30 %, at sub-degree 
resolution and the large-angular scales of both the CMB 
temperature anisotropics and the temperature and polarization 
diffuse emission from Galactic dust. Maps of the tempera- 
ture diffuse Galactic dust emission are available at the four 
Archeops frequency channels, 143, 217, 353 and 545 GHz. 
These are very useful maps as intermediate resolution products 
between the FIRAS and the expected Planck HFI maps 
l |Bemard (2004)| l. Foreground-cleaned CMB maps were pro- 
duced for the lowest frequency channels at 143 and 217 GHz 
using the information on dust and atmospheric emission 
provided by the high frequency ones. These maps provided 
the first simultaneous determination of the Sachs-Wolfe 
plateau and of the first acoustic peak of the CMB anisotropics 
temperature power spectrum ( Benoit e t al. (200 3a)[l and more 
recently of the second acoustic peak jTristari^^l. (200 5b)). 
By combining those maps with the data from the WMAP 
experiment and the 2MASS catalog of galaxies we ob- 
tained a local statistical detection of the SZ effect in clusters 
CHernandez-Monteagudo et al. (2006) ). In a forthcoming paper 
(Patanchon et al. (2005)| ), an analysis of the level of any diffuse 
SZ emission will be presented, using the Archeops and WMAP 
data in order to have a broad electromagnetic spectral leverage. 
Finally, /, Q and U maps of the polarized diffuse emission 
of Galactic dust were constructed at 353 GHz combining the 
measurements from the six polarized sensitive bolometers. 
These maps allowed us, for the first time, to characterize the 
polarized diffuse emission from Galactic dust in the Galactic 
plane (Benoit et al. (200 4)^ and to estimate the polarized power 
spectra of the diffuse Galactic dust emission at intermediate 
and high Galactic latitudes |Ponthieu et al. 92005)| ). 

The Archeops data are the first available data which present 
very similar characteristics to those of the Planck HFI instru- 
ment. This is because the instrumental configuration, the acqui- 
sition system and the scanning strategy in Archeops and Planck 
are very similar. There are few important differences between 
these two data sets as for example the presence of an atmo- 
spheric signal in the Archeops data which would not be at all 
present in the Planck data and which is one of the most impor- 
tant systematics in Archeops. However, in many other aspects 
they are sufficiently similar to consider that the techniques and 
methods applied for the processing of the Archeops data will 



be of great use for processing the Planck HFI data. In this 
sense, the processing of the Archeops data was for us a learning 
process towards the analysis of the Planck HFI data. Actually, 
most of the preprocessing, decorrelation, deconvolution from 
the bolometer time constant, beam pattern reconstruction, noise 
spectrum estimation, destriping, calibration and power spec- 
trum estimation methods are currently being adapted to the 
Planck data within the Planck HFI Level 2 data processing. 

We encourage interested parties to contact members of the 
Archeops collaboration to any specific scientific project (for 
example correlation with other data sets) using the Archeops 
data. 

Appendix A: Intercalibration procedure 

We describe here the details of the intercalibration procedure 
of the Archeops bolometers discussed in Sect. 19.41 

A.0.4. Modeling 

Let's formulate the problem as follows : s\{b), . . . , SN{b) are 
profiles (e.g. of the Galaxy), measured by different bolome- 
ters; b stands for the Galactic latitude, and runs from 1 to B\ Sb 
is the estimated profile, which has to be determined. 
We can write the estimated profile as : 

s jib) = a j-^(b) + tijib) (A.l) 

where nj{b) is the noise in the bin b of profile j and aj the 
associated calibration coefficient. 

If we assume Gaussian white noise, we can write the as 

,2 V V (sj(b) - aj^b))^ 
^ "hi hi TT^i ^^^^ 

with o-j(b)^ the noise variance {nj(b)^) (we neglect noise 
correlations between detectors). 

A.0.5. Constraint 

We notice that the is invariant under the following transfor- 
mation : 

Skib) ^ /5-Mb) 

This degree of freedom is due to the fact that we can only 
determine the intercalibration coefficient up to a constant fac- 
tor. We must choose a constraint on the parameters in order to 
converge to a unique solution of the equations. 

Let's choose a general relation : g{{aj}, [sib)}) - 0. Using 
the method of Lagrange's multiplier, we thus have to minimize 
the function 
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+ Ag{{aiU-s(b)]) 



(A.4) 



with respect to {ff,), and A. The conditions of mini- 
mum lead to the three equations : 



^ = 0^g({a,),{^W)) = 

OA 



dx_ 



aiisiib) - ais(b)) 



+ A- 



ds(b) 



= 



dai ^ 



7s(b)(si(b) - aisjb)) dg 
erf, dUi 



(A.5) 
(A.6) 

(A.7) 



Multiplying Eq. ( IA.6> by lib) and summing over b, and 
multiplying Eq. ( IA.7> and summing over /, we find that A - Q 
is the only solution. We thus find the two following relations : 



-s(b) 



a-j(b)^ 

2; 



and 



y s(b).ik(h) 



whatever is the constraint. 

These equations can be solved by iteration : starting from 
any set of {a,), we can calculate lib) with Eq. ( IA.8> . and cal- 
culate new {a,} with Eq. (IA.9> . At each step of the iteration, 
we have to check that the constraint is satisfied, or to impose 
it. The iteration ends when the relative differences between two 
successive steps is small enough. 

A.0.6. Error matrix 

We can develop any functions around its minimum as : 



Since F is a positive definite matrix, it can be diagonalized, 
with all its eigen values positive. Let's call M the change of 
frame matrix, so that F - M'DM, where D is diagonal, and 
A0' = MA0. M is orthogonal, so that M"' = M'. The correla- 
tion between two measurements will be given by : (A0A0') - 
(M'AQ'A&'M) = M'<A0'A0">M. The central part is the cor- 
relation matrix in the diagonal frame : it can be calculated di- 
rectly to give (A0'A0") - D^^. We then deduce the corre- 
lation matrix in the original frame : <A0A0') = M'D'^M - 
(M'DMY^ - All this calculation is made with the as- 
sumption that the likelihood is very close to a Gaussian around 
its maximum. We will see in the following that it is the case for 
our particular case. 

In our particular case, the Fisher matrix F is dimension (A^-H 
B-i-l)x(A^H-B-i-l), and we can compute it analytically for any 
constraint g : 



1 d^x^ 



2 dA^ 
1 5V 



= 



dg 



2 dAdl{b) ds(b) 

1 

2 dAdoj dai 



(A.8) L^JL- = 2L 



BY 



2 dl{b)dl{q) 



(A.9) 1 5V 



2 ds{b)dak 



1 d'x' 



2aks(b) - Sk(b) 

CTkib)^ 

~sibf , 



a-iib)^"'^ 



(A. 12) 
(A.13) 
(A. 14) 

(A. 15) 

(A. 16) 
(A. 17) 



2 doida i , , 

■' b-l 

(taking all the parameters at the minimum, including A 

0). 

Fig. A.l. Simulation profiles 



xH&) = x'i^^"""^) 

2 50,5©/ ' ' 'y^' j > 
+ 0(&^) (A. 10) 

where is the vector of parameters and 0("""' is the mini- 
mum of the;^'^. The matrix j]^^ is called the Fisher matrix 
(noted F in the following) and is the inverse of the matrix of 
correlation of the parameters, as can be easily seen : calling A0 
the vector 0-0*"""', the probability that the real parameters are 
can be written using the likelihood function X = exp(-;^f^/2), 
i.e. : 

£oce —. (A.ll) 



A. 0.7. Robustness tests 

We have performed a bunch of tests in order to validate this 
intercalibration method. First, we have tested by Monte Carlo 
simulations the reliability of the error bars computed using the 
Fisher matrix. Second, we have compared the influence of the 
choice of the constraint on the intercalibration coefficients and 
their errors : when comparing what is comparable, i.e. the ratio 
of intercalibration coefficients, we found no difference neither 
in the value nor in the error. We have chosen to impose the 
constraint g{{'s(b),{aj}}) - ai - I (i.e. the first profile has a 
relative calibration coefficient with respect to the average pro- 
file of 1). Finally, we have compared the iterative minimization 
method (using alternatively Eas. lA.8l and lA.9t with a standard 
minimization program (Minuit, from the CernLib). Differences 
were below the numerical precision. 
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Fig. 1. Adiabatic exponent Fi. Fi is plotted as a function of Ig internal energy [ergg '] and Ig density [g cm ^]. 



Introduction 

In the nucleated instability (also called core instability) hypothesis of giant planet formation, a critical mass for static core 
envelope protoplanets has been found. Mizuno ( 1980 1 determined the critical mass of the core to be about 12 (M® = 5.975 x 
10^^ g is the Earth mass), which is independent of the outer boundary conditions and therefore independent of the location in the 
solar nebula. This critical value for the core mass corresponds closely to the cores of today's giant planets. 
■ ■ Although no hydrodynamical study has been available many workers conjectured that a collapse or rapid contraction will 
, ensue after accumulating the critical mass. The main motivation for this article is to investigate the stability of the static envelope 
' at the critical mass. With this aim the local, linear stability of static radiative gas spheres is investigated on the basis of Baker's 
\ ( I1966t standard one-zone model. 

' Phenomena similar to the ones described above for giant planet formation have been found in hydrodynamical models con- 
. ^ \ cerning star formation where protostellar cores explode (Tscharnuter 1 1 9 87l Balluch[r988 1, whereas earUer studies found quasi- 
steady collapse flows. The similarities in the (micro)physics, i.e., constitutive relations of protostellar cores and protogiant planets 
serve as a further motivation for this study. 



(N 



1 . Baker's standard one-zone model 



In this section the one-zone model of Baker J1966> . originally used to study the Cepheid pulsation mechanism, will be briefly 
\^ reviewed. The resulting stability criteria will be rewritten in terms of local state variables, local timescales and constitutive 
\^ relations. 

Baker (I1966> investigates the stability of thin layers in self-gravitating, spherical gas clouds with the following properties: 

o ; 

. - hydrostatic equilibrium, 
\ - thermal equilibrium, 
f~i . - energy transport by grey radiation diffiision. 

^: 

' For the one-zone-model Baker obtains necessary conditions for dynamical, secular and vibrational (or pulsational) stability (Eqs. 
^ (34a, b, c) in Baker 1966) . Using Baker's notation: 

Mr mass internal to the radius r 
I ' m mass of the zone 



ro unperturbed zone radius 
po unperturbed density in the zone 



To unperturbed temperature in the zone 
■ - - ' L,-o unperturbed luminosity 

E(h thermal energy of the zone 

and with the definitions of the local cooling time (see Fig.QJ 

^th 

Tco = — , (1) 

and the local free-fall time 



Tff 



— ^ (2) 

V 32G 3Mr ' ^ ^ 



Baker's K and ctq have the following form: 

C.0 = ^- (3) 
V8 ^ff 

/.= ^ii; (4) 

TT Tea 

where Eih ~ m(f o/Po) has been used and 
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Table 1. Opacity sources. 



Source 


Tim 


Yorke 1979, Yorke 1980a 


< 1700" 


Krugel 1971 


1700 <T < 5000 


Cox & Stewart 1969 


5000 < 



" This is footnote a 



is a thermodynamical quantity which is of order 1 and equal to 1 for nonreacting mixtures of classical perfect gases. The physical 
meaning of ctq and K is clearly visible in the equations above, ctq represents a frequency of the order one per free-fall time. K 
is proportional to the ratio of the free-fall time and the cooling time. Substituting into Baker's criteria, using thermodynamic 
identities and definitions of thermodynamic quantities. 



{ d\nP \ 



X p 



d\nP\ 



d\nP 



'ad 



f^lnT 
[d\nP 



d\nP 
d\nT 



d\nK 
d\nT 



one obtains, after some pages of algebra, the conditions for stability given below: 



^4(3ri -4) > 

8 tI 



-r,v 



3 

4 TcoT^- 



1 Vad 



1 - 3/4;rp 



Xt 



{Kj -A) + Kp+\ 



4Vad - (Vad/Cr + Kp) - — 

31 1 



> 

> 



(6) 
(7) 
(8) 



For a physical discussion of the stability criteria see Baker (I1966> or Cox (I1980> . 

We observe that these criteria for dynamical, secular and vibrational stability, respectively, can be factorized into 

1 . a factor containing local timescales only, 

2. a factor containing only constitutive relations and their derivatives. 

The first factors, depending on only timescales, are positive by definition. The signs of the left hand sides of the inequalities (|6}, 
Q and therefore depend exclusively on the second factors containing the constitutive relations. Since they depend only on 
state variables, the stability criteria themselves are functions of the thermodynamic state in the local zone. The one-zone stability 
can therefore be determined from a simple equation of state, given for example, as a function of density and temperature. 
Once the microphysics, i.e. the thermodynamics and opacities (see Table[0, are specified (in practice by specifying a chemical 
composition) the one-zone stability can be inferred if the thermodynamic state is specified. The zone - or in other words the layer 
- will be stable or unstable in whatever object it is imbedded as long as it satisfies the one-zone-model assumptions. Only the 
specific growth rates (depending upon the time scales) will be diff'erent for layers in different objects. 

We will now write down the sign (and therefore stability) determining parts of the left-hand sides of the inequalities (|6}, Q 
and and thereby obtain stability equations of state. 

The sign determining part of inequality (|6j is 3ri - 4 and it reduces to the criterion for dynamical stability 



(9) 



Stability of the thermodynamical equilibrium demands 

Xp > 0, c, > 0, 

and 



(10) 



(11) 
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Fig. 2. Vibrational stability equation of state 5'vib(lg e, Igp). > means vibrational stability. 



holds for a wide range of physical situations. With 

r3 - 1 = 4— > ^12) 

pT Ci, 

ri^Xp+Xji^-i) > (13) 
r3 - 1 

V,d = > (14) 

i 1 

we find the sign determining terms in inequalities and (jSjl respectively and obtain the following form of the criteria for 
dynamical, secular and vibrational stability, respectively: 

3ri-4=:5dy„> (15) 

1 - 3/4;rp 

-(Kj-4) + Kp+l-.S,^c> (16) 

4 

4Vad - (Vad/fr +Kp)-— =: 5vib > . (17) 
il I 

The constitutive relations are to be evaluated for the unperturbed thermodynamic state (say (po, Tq)) of the zone. We see that 
the one-zone stability of the layer depends only on the constitutive relations Fi, Vad, Xt^ Xp^ i^p' k^t- These depend only on the 
unperturbed thermodynamical state of the layer. Therefore the above relations define the one-zone-stability equations of state 
Sdyn, Ssec and 5vib- See Fig.|2lfor a picture of Snb- Regions of secular instability are listed in Table 1. 



2. Conclusions 

1. The conditions for the stability of static, radiative layers in gas spheres, as described by Baker's (19661 standard one-zone 
model, can be expressed as stability equations of state. These stabihty equations of state depend only on the local thermody- 
namic state of the layer. 

2. If the constitutive relations - equations of state and Rosseland mean opacities - are specified, the stability equations of state 
can be evaluated without specifying properties of the layer. 

3. For solar composition gas the A:-mechanism is working in the regions of the ice and dust features in the opacities, the H2 
dissociation and the combined H, first He ionization zone, as indicated by vibrational instability. These regions of instability 
are much larger in extent and degree of instability than the second He ionization zone that drives the Cepheid pulsations. 
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